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PREFACE 


This  final  technical  report  describes  the  results  of  one  and  a  half  years  of  effort  on  the 
Integrated  Optical  Transceiver  program  to  develop  100GHz  operation  in  electro-optic 
transceivers.  The  program  resulted  in  the  development  of  new  concepts  in  device 
operation  and  fabrication  and  growth  for  a  unified  set  of  components  in  the  Inversion 
Channel  Technology.  The  designs  analyses,  fabrication  sequence  design  and  photomask 
design  are  described  in  detail. 


1.  SUMMARY 


The  integrated  Optical  transceiver  program  was  initiated  on  Sept.  16, 1995  and 
was  concluded  on  March  31,  1997.  This  phase  I  effort  was  concerned  with  establishing 
the  feasibility  of  implementing  an  electro-optic  data  link  with  a  digital  data  capability  of 
lOOGb/s  or  an  analog  transmission  rate  of  100GHz  .  To  achieve  this  goal,  three  new 
concepts  were  introduced  which  included  a  grating  coupled  vertical  cavity  laser  and 
detector,  a  traveling  wave  electro-optic  device  design  that  builds  upon  the  existing 
inversion  channel  technology,  and  an  implementation  of  the  inversion  channel  bipolar 
transistor  that  can  achieve  very  low  parasitic  resistance  by  virtue  of  its  inverted  structure. 
The  effort  in  this  program  established  a  technology  plan  to  achieve  100GHz  operation  in 
traveling  wave  operation.  The  principles  of  diffraction  in  second  order  corrugated 
waveguides  were  established  and  tools  were  developed  to  design  the  grating  coupled 
laser.  The  principles  of  velocity  matching  in  electro-optic  devices  was  developed  through 
the  analysis  of  the  grating  coupled  detector  for  high  speed  operation.  The  fabrication 
sequence  of  the  Inversion  Channel  Technology  was  completely  revised  to  enable 
waveguides  and  traveling  wave  operation.  Growth  and  fabrication  of  these  concepts  is 
the  next  step. 
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2.  BACKGROUND 


The  integrated  optical  transceiver  program  was  designed  to  determine  the 
feasibility  of  a  high  bandwidth,  frequency-modulated  transceiver  for  optical  data 
communications  to  and  from  an  electronic  source.  The  E/O  and  O/E  signal  conversions 
need  to  be  suitable  for  digital  transmission  at  100  Gb/s  with  a  bit  error  rate  of  10'12  or 
analog  transmission  with  a  3  dB  cut-off  frequency  of  100  GHz.  Our  interpretation  of  the 
program  requirements  was  that  a  single  channel  operating  at  100  GHz  was  desired  and 
therefore  the  goal  was  operation  of  electronic  amplifying  devices  with  cut-off  frequencies 
of  150  GHz-200  GHz  and  optoelectronic  devices  with  cut-off  frequencies  of  100  GHz. 
The  higher  frequency  operation  of  transistors  is  needed  to  build  circuits  with  this 
performance.  Multiple  channels  could  be  considered  to  reduce  the  bandwidth  requirement 
on  the  emitters  and  detectors  but  it  would  still  require  the  same  speed  of  electronics. 

An  analysis  of  these  requirements  indicates  that  the  problems  of  interconnection 
of  the  electronic  and  optoelectronic  components  are  major  obstacles  to  these  goals  and 
that  a  monolithic  approach  that  integrates  all  of  these  components  on  the  same  substrate  is 
the  only  practical  way  to  approach  this  performance.  Above  20  GHz  the  problems  of  wire 
or  even  bump  bonding  of  transistor  circuits  to  opto  devices  becomes  very  difficult.  The 
metal  traces  on  the  chips  must  be  matched  to  50  ohms  and  the  solder  bump  impedance 
must  appear  as  50ohms.  If  wire  bonding  and  MM  wave  connectors  are  used  these  must 
also  be  matched  to  the  metal  traces.  To  achieve  a  manufacturable  (affordable)  result 
monolithic  integration  is  the  only  solution. 

Our  proposal  to  achieve  these  goals  was  based  upon  a  monolithic  integrated 
technology  approach  which  was  developed  by  WPAFB  through  funding  to  AT&T  over 
the  period  from  1987-1994.  This  technology  has  been  designated  as  the  Inversion 
Channel  Technology  (ICT)  because  it  utilizes  the  inversion  channel  formed  at  a 
modulation  doped  interface  as  the  active  element  in  transistor,  laser,  detector,  and 
modulator  devices.  It  also  produces  all  of  these  devices  with  a  single  unified  fabrication 
sequence  which  uses  the  same  set  of  processing  steps  to  realize  all  components 
simultaneously.  The  starting  wafer  is  grown  by  MBE  and  a  single  layer  sequence  is  used 
for  all  devices  without  the  requirement  for  any  regrowth  steps.  This  technology  is  actually 
an  optoelectronic  version  of  the  very  successful  PHEMT  technology  because  it  employs 
essentially  the  same  growth  structure.  Both  wafer  growths  utilize  a  modulation  doped  III- 
V  interface  in  which  a  multiple  quantum  well  of  strained  InGaAs  forms  a  channel  at  a 
heterostructure  interface.  Both  growths  start  from  a  SI  GaAs  substrate  and  place 
moderately  p-doped  AlGaAs  below  the  quantum  wells.  In  these  respects  they  are 
essentially  identical.  Where  they  differ  is  above  the  quantum  well(s).  Whereas  the 
PHEMT  grows  only  a  thin  layer  of  undoped  AlGaAs  upon  which  a  schottky  gate  is 
formed,  the  opto  PHEMT  grows  layers  of  p  doped  AlGaAs  which  are  terminated  in  p+ 
doping  to  form  an  ohmic  contact.  It  is  this  innovation  which  allows  the  formation  of  the 
laser,  the  detector  and  the  transistor  simultaneously. 

The  ICT  technology  was  developed  at  AT&T,  first  with  edge  emitting  optical 
devices  and  then  with  vertical  cavity  optical  devices.  The  vertical  cavity  technology 
version  was  barely  to  the  point  of  full  demonstration  when  this  program  was  interrupted 
by  AT&T  and  the  PI  moved  to  the  University  of  Connecticut.  What  had  been 
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demonstrated  was  the  operation  of  the  transistor,  the  switching  laser  and  the  detector  all 
from  the  same  epitaxial  wafer.  In  terms  of  performance  at  that  point  in  time,  the  transistor 
operated  at  14  GHz  using  Be  for  p  doping  for  a  1pm  device  gate  length  but  when  Carbon 
was  used  for  the  p  doping  in  the  growth  and  e-beam  lithography  was  used  to  generate  0.5 
pm  gates,  we  obtained  40  GHz  cut-off  frequencies.  The  vertical  cavity  laser  3dB 
frequency  was  3.8  GHz.  The  resonant  cavity  detector  from  this  epi  growth  was 
demonstrated  with  near  unity  quantum  efficiency  but  speed  measurements  were  not 
obtained.  These  component  demonstrations  were  made  from  the  same  epitaxy  and  using 
the  same  fabrication  sequence.  Due  to  the  success  of  the  individual  experiments, 
photomasks  were  designed  to  produce  functional  opto  IC’s.  Initial  fabrication  had  begun 
and  a  couple  of  chips  made  it  almost  to  the  end  of  the  sequence.  These  chips  however  did 
not  qualify  for  deposition  of  dielectric  stacks  to  demonstrate  the  opto  devices  and  so  final 
testing  of  full  IC’s  had  only  just  begun.  Testing  was  restricted  to  the  FET  electronic 
devices  where  divide  by  two  circuits  and  flip  flops  were  demonstrated  and  to  individual 
laser  devices  where  a  threshold  of  2  mA  for  the  DOES  laser  with  proton  implant  isolation 
was  achieved. 

We  have  proposed  for  the  integrated  optical  transceiver  program  to  use  the  ICT  as 
a  starting  point  to  realize  monolithically  integrated  components.  However,  to  have  a  hope 
of  achieving  the  very  high  bandwidth  operation,  several  device  innovations  were  required 
because  there  are  no  known  device  structures  either  within  the  ICT  or  elsewhere,  that 
could  achieve  100GHz  operation  in  the  laser.  There  are  also  no  known  detector  structures 
that  could  achieve  these  bandwidths  within  the  constraints  of  a  monolithic  integrated 
circuit.  In  addition,  although  the  laser  and  detector  were  not  required  to  be  part  of  the 
same  chip,  our  concept  was  to  produce  one  integrated  circuit  which  would  perform  the 
transceiver  function  at  either  end  of  an  optical  link  so  the  integration  would  include  the 
laser,  detector  and  FET  electronics  .  There  were  therefore,  several  challenges  to  be  met  in 
our  program  and  these  are  listed 

•Devise  a  laser  structure  within  the  ICT  to  achieve  1 00  GHz  operation.  This 
means  reworking  the  DOES  laser  or  HFET  laser  concept  to  circumvent  the  conventional 
limits 

•Devise  a  detector  within  the  ICT  to  achieve  100  GHz  operation.  The  HFET 
detector  concept  must  also  be  reworked  to  address  the  high  speed  requirement. 

•The  new  concepts  required  on-chip  waveguides.  Therefore  the  processing  flow 
sequence  to  build  the  ICT  had  to  be  totally  redesigned.  As  it  turned  out ,  the  photomask 
sequence  and  usage  were  altered  substantially  from  the  final  AT&T  version  in  order  to 
accommodate  the  new  concepts.  Therefore  a  major  mask  redesign  had  to  be  implemented 

•The  new  laser  and  detector  designs  required  CAD  design  tools  for  gratings  that 
did  not  exist  and  thus  a  major  effort  was  needed  to  analyze  and  understand  the 
performance  parameters  of  second  order  gratings. 

•To  build  these  complex  structures  requires  a  state-of-the-art  fabrication  and 
growth  facility.  A  major  burden  on  the  PI  since  arriving  at  UConn  was  to  put  this 
capability  in  place.  The  task  was  finished  only  as  the  program  drew  to  a  close  and  was  the 
major  reason  why  so  little  progress  was  made  within  the  one  year  on  the  reduction  to 
practice  of  these  concepts. 
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In  the  next  sections  we  describe  the  innovations  ,  simulations ,  analyses  and 
designs  that  were  implemented  resulting  in  the  final  photomasks  that  were  made  for  the 
demonstration.  We  will  also  describe  our  efforts  to  reach  a  demonstration  and  how 
funding  limitations  have  curtailed  our  attempts  to  move  forward.  Our  present  efforts  to 
acquire  these  funds  through  other  channels  and  the  current  status  and  plans  will  be 
mentioned  briefly. 
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3.  LASER  DESIGN  REQUIREMENTS  FOR  100GHz  OPERATION 


The  laser  structure  to  achieve  100  GHz  operation  must  be  a  vertical  cavity  device 
because  it  is  only  the  vertical  cavity  in  which 

(1)  the  separate  confinement  heterostructure  (SCH)  region  can  be  made  very  short 
because  optical  confinement  by  the  SCH  region  is  not  required  (the  optical 
confinement  in  the  vertical  cavity  is  the  ratio  of  the  quantum  wells  to  the  cavity 
length).  When  the  SCH  region  is  very  short  (a  few  hundred  angstroms),  the  low 
frequency  roll-off  effects  are  minimized  so  that  low  frequency  roll-off  can  be 
extended  to  beyond  100GHz. 

(2)  The  photon  lifetime  of  the  cavity  can  be  made  veiy  short.  Since  the  stimulated 
lifetime  tracks  the  photon  lifetime  to  within  1/Q  (Q  is  the  quality  factor  of  the  laser) 
then  this  determines  the  ultimate  laser  speed.  In  calculations  of  the  response  the 
parameter  of  interest  is  the  differential  stimulated  lifetime.  In  conventional  laser 
jargon  this  is  proportional  to  the  differential  gain. 

When  both  these  effects  have  been  addressed  for  high  speed  operation  then  the 
remaining  limitations  on  the  laser  response  are  the  time  constants  presented  by  resistance 
and  capacitance.  First  there  is  the  time  constant  determined  by  the  resistance  and 
capacitance  internal  to  the  laser  diode.  In  the  vertical  cavity  laser,  the  resistance  is 
primarily  the  contact  resistance  of  the  p  type  contact  and  its  associated  current  crowding 
and  the  capacitance  is  that  of  the  i  region  of  the  pin  structure.  The  capacitance  is  fixed  by 
design  but  the  resistance  can  be  minimized  by  innovations  in  the  metallurgy  and  the 
geometry  of  the  contact  areas.  Reducing  this  resistance  becomes  the  ultimate  limitation  in 
increasing  the  bandwidth  of  the  laser.  The  second  time  constant  is  the  one  determined  by 
the  laser  capacitance  and  the  load  impedance  driven  by  the  laser.  For  discrete  laser 
devices  this  is  always  50  ohms.  In  order  to  address  this  issue,  a  traveling  wave  type  of 
operation  is  required.  A  concept  was  therefore  introduced  by  which  the  virtues  of  the 
vertical  cavity  laser  could  be  maintained  within  the  context  of  a  traveling  wave  device. 
The  physical  implementation  of  the  concept  was  within  the  ICT  and  the  device  is  shown 
in  cross-section  along  the  direction  of  the  traveling  wave  in  Fig.l  and  in  a  perspective 
view  from  the  input  port  in  Fig.2. 
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Fig.  1  Grating  Coupled  Traveling  Wave  Laser 


U2  4/2 


0.98  m  In  s.l.  GaAs 


Fig.2  Traveling  wave  optoelectronic  device  cross-section  in  the  ICT 
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4.  THE  GRATING  COUPLED  VERTICAL  CAVITY  LASER 


The  innovation  is  to  introduce  a  diffraction  grating  within  the  vertical  cavity  as  shown  in 
Fig.  1 .  The  grating  is  formed  in  the  first  layer  of  the  mirror  stack.  It  is  a  blazed  grating 
which  diffracts  light  preferentially  in  one  direction.  The  DBR  mirrors  of  the  vertical 
cavity  device  are  deposited  on  top  of  the  mirror  and  of  course  the  bottom  mirror  is  formed 
under  the  device  during  the  growth.  We  call  this  device  the  Grating  Coupled  Vertical 
Cavity  Surface  Emitting  Laser  (GC  VCSEL).  To  understand  the  operation  of  the  GC 
VCSEL  consider  the  VCSEL  operation  at  the  back  end  of  the  device  (x  =  0  in  Fig.l) .  The 
front  end  at  x  =  L  is  the  output  port  where  light  is  to  be  emitted.  The  current  is  to  be 
applied  uniformly  along  the  length  of  the  device.  When  current  is  applied,  light  is 
generated  vertically  in  the  VCSEL  at  x  =  0.  Consider  the  optical  wave  in  the  cavity  as  it 
passes  through  the  grating.  A  small  fraction  is  diffracted  to  the  left.  Due  to  the  asymmetry 
of  the  blaze  a  much  smaller  fraction  is  diffracted  to  the  right  (this  is  illustrated  by  the 
relative  sizes  of  the  arrows  in  Fig.l).  Now  if  the  grating  were  of  zero  thickness,  then  the 
remainder  of  the  light ,  which  is  most  of  the  energy,  would  pass  through  the  grating,  be 
reflected  back  and  pass  through  the  grating  in  the  reverse  direction.  Due  to  the 
reversibility  of  the  diffraction  process,  the  majority  of  this  light  would  be  diffracted  in  the 
reverse  direction  with  the  net  result  that  the  light  would  be  coupled  fairly  equally  in  both 
directions.  However,  the  grating  has  a  finite  thickness  and  in  addition  the  penetration  of 
the  light  into  the  mirror  is  small  for  a  high  reflectivity  mirror.  Thus  the  intensity  of  the 
light  traveling  downward  into  the  laser  after  reflection  from  the  mirror  decreases 
dramatically  with  penetration  into  the  mirror.  It  follows  then  that  the  amplitude  of  the 
light  incident  on  the  grating  from  within  the  cavity  is  substantially  greater  than  the 
amplitude  of  the  light  incident  on  the  grating  by  reflection  from  the  mirror.  Therefore  the 
fraction  of  light  diffracted  to  the  right  is  substantially  less  than  the  fraction  of  light 
diffracted  to  the  left. 

The  design  criterion  that  is  followed  is  that  the  thickness  of  the  diffraction  grating 
should  be  approximately  equal  to  the  penetration  depth  of  the  mirror  ( this  is  also 
determined  by  the  standing  wave  effects  in  the  grating  which  are  found  from  the 
diffraction  analysis).  This  situation  is  the  one  that  will  be  implemented  in  practice 
because  for  a  high  reflectivity  mirror,  the  decay  of  intensity  to  about  10%  occurs  within 
the  first  14  wavelength  of  the  mirror  stack  and  additionally,  it  is  most  practical  to  form 
the  grating  by  etching  through  approximately  the  first  layer  of  the  mirror  stack.  In  the 
fabrication  sequence,  the  dielectrics  chosen  to  form  the  mirrors  are  SiC>2  and  undoped 
sputtered  GaAs.  This  choice  is  dictated  by  the  very  large  index  difference  which  reduces 
the  number  of  required  pairs  and  the  fact  that  our  mirror  must  be  deposited  during  the 
device  processing.  Given  these  layer  components  we  have  two  choices  to  implement  the 
grating  l)deposit  a  layer  of  GaAs  and  then  pattern  and  etch  this  layer  followed  by  a  14 
wave  of  SiC>2  and  then  GaAs  etc.  2)deposit  a  14  wave  layer  of  SiC>2,  pattern  and  etch  this 
layer  and  then  deposit  a  14  wave  of  GaAs  and  then  SiCh  and  then  GaAs  etc.  The  second  of 
these  is  the  one  we  will  start  with  in  the  fabrication.  The  choice  of  14  wave  thickness  is 
also  supported  by  the  grating  efficiency  and  its  dependence  on  the  index  difference 
between  the  layers  of  the  stack  as  discussed  in  the  next  section  on  grating  efficiencies.  In 
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Fig.l,  the  relative  sizes  of  the  arrows  indicate  that  the  majority  of  light  is  diffracted  to  the 
right. 

5.  DETERMINATION  OF  DIFFRACTION  EFFICIENCIES  IN  THE  GRATING 
COUPLED  VERTICAL  CAVITY  LASER 

The  quantitative  prediction  of  diffracted  power  has  been  a  subject  of  much  study 
over  the  years  primarily  because  of  its  importance  to  the  operation  of  the  DFB  laser.  In 
these  cases  the  gratings  are  usually  designed  to  be  first  order  since  it  is  the  first  order 
forward  and  backward  traveling  waves  which  form  the  basis  of  coupled  mode  theory. 
There  are  some  instances  where  second  order  has  been  used  in  the  DFB  to  produce  an 
optical  loss  mechanism  in  an  effort  to  stabilize  the  mode  position  as  described  by 
Kazarinov  and  Henry  [6](however,  the  laser  output  is  still  via  the  first  order  wave).  More 
generally  the  second  order  grating  has  been  used  as  the  output  reflector  for  DBR  lasers  in 
order  to  obtain  vertical  emission.  Since  the  second  order  grating  when  implemented  as  a 
waveguide  corrugation  produces  both  first  order  diffraction  in  the  guide  and  second  order 
diffraction  normal  to  the  guide  then  it  can  be  used  in  the  DBR  to  provide  both  reflection 
for  the  guided  wave  and  output  coupling  of  the  laser  light. 

At  the  beginning  of  this  work,  the  best  available  description  of  waveguide 
diffraction  was  the  work  of  Streifer  [5]  and  Kazarinov  [6].  All  of  the  existing  work 
describes  the  diffraction  in  a  general  way  which  includes  all  possible  diffraction  orders. 
The  results  are  complicated,  difficult  to  understand  and  require  a  full  numerical  analysis 
to  implement  to  obtain  the  quantitative  predictions  required  for  this  work.  We  needed  a 
simpler  design  tool  with  the  same  level  of  accuracy.  Therefore  it  was  undertaken  to  solve 
this  problem  in  a  more  direct  manner  to  produce  more  useable  results.  We  are  only 
interested  in  ideal  second  order  diffraction  so  this  was  used  as  a  starting  point.  It  turns  out 
that  for  this  case  the  equations  can  be  simplified  and  only  one  differential  equation  needs 
to  be  solved  for  a  symmetrical  grating.  With  this  approach  a  simple  predictor  for 
diffraction  efficiency  was  produced  for  both  exact  or  approximate  second  order  gratings, 
i.e.  the  approximate  second  order  grating  is  one  in  which  the  grating  pitch  is  either 
slightly  smaller  or  larger  than  the  exact  value.  This  solver  runs  on  a  200  MHz  PC  in  about 
150  secs  to  produce  a  complete  calculation  of  efficiency  versus  grating  depth  for 
example.  The  basis  for  the  model  was  a  three  layer  waveguide  characterized  by  indices 
nj,  n2,  and  n3  in  the  regions  above  the  guide,  in  the  guide  and  below  the  guide.  This 
analysis  was  performed  for  diffraction  from  the  guide  to  the  direction  normal  to  the  guide. 
This  work  was  published  and  a  copy  of  the  publication  “Determination  of  Diffraction 
Efficiency  for  a  Second-Order  Corrugated  Waveguide”  is  included  in  appendix  A. 

All  of  the  published  work  on  waveguide  diffraction  has  been  for  the  process  of 
diffraction  of  the  guided  wave  into  an  unguided  (leaky)  wave  propagating  away  from  the 
guide.  Of  most  interest  from  the  point  of  view  of  the  grating  coupled  VCSEL,  is  the 
diffraction  of  light  from  a  normally  propagating  wave  into  a  guided  wave  (here  the  wave 
is  generated  within  the  vertical  cavity  and  propagates  into  the  grating  as  shown  in  Fig.3  ). 
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Fig.3  Topologies  for  the  diffraction  of  light  normal  to  a  corrugated  waveguide: 
left,  light  incident  externally  upon  a  corrugated  waveguide 
right,  light  incident  upon  corrugation  from  within  guide  (w  or  w/o  grating) 
This  is  the  reverse  to  the  process  described  in  the  above  paper  but  the  efficiency  of  the 
process  is  not  necessarily  the  same  in  spite  of  the  fact  that  optical  propagation  of  energy 
through  a  passive  structure  should  be  reversible.  This  is  because  of  the  transmission  of 
energy  straight  through  the  guide  which  is  not  present  in  the  out-of-guide  case.  Therefore 
this  mechanism  was  also  characterized  using  the  same  basic  mathematical  approach.  This 
calculation  was  also  performed  for  the  generic  three  layer  waveguide  with  a  blazed 
grating.  This  work  is  summarized  in  the  paper  “Diffraction  Into  a  Corrugated  Waveguide 
From  Normally  Incident  Radiation”  which  has  been  submitted  for  publication  and  forms 
appendix  B.  One  further  piece  of  necessary  work  was  the  determination  of  the  optimum 
blaze  angle  to  maximize  the  diffraction  in  one  direction  and  to  minimize  it  in  the  reverse. 
This  design  tool  has  also  been  developed  and  the  results  will  be  submitted  for  publication 
soon.  Some  typical  results  are  shown  in  Fig.4  where  we  plot  the  diffraction  efficiency  into 
the  guide  for  a  corrugated  waveguide  as  the  blaze  angle  is  varied.  It  shows  the  optimum 
angle  for  typical  index  parameters(  ni=n3=3.4  and  n2=3.6)  is  about  41°  and  the  ratio  of  the 
diffraction  in  the  two  directions  is  about  135. 


Fig.  4  Diffraction  efficiency  versus  grating  angle  in  a  three-layer  dielectric  guide 
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The  above  two  articles  have  produced  the  tools  needed  for  the  design  of  the 
GCVSEL.  However ,  in  the  GCVSEL  the  waveguide  cladding  is  formed  by  the  Bragg 
mirrors  and  thus  the  structure  forms  a  multilayer  waveguide  with  very  complicated 
propagation  coefficients.  In  order  to  apply  the  diffraction  model,  it  is  necessary  to 
represent  the  VCSEL  waveguide  by  an  equivalent  three-layer  guide.  To  do  this  it  is 
necessary  to  represent  the  top(bottom)  stack  and  its  interface  to  air  (substrate)  by  an 
infinite  region  with  a  single  index.  Since  the  material  between  the  mirrors  can  be 
represented  by  a  single  index  determined  as  the  average  over  its  layers,  then  the  result 
will  be  a  three-layer  waveguide.  The  wave  propagating  in  a  vertical  cavity  laser  is 
characterized  by  a  penetration  depth  .  This  is  normally  used  to  determine  the  effective 
cavity  length  of  the  VCSEL  as 

^eff  =  +  ^pb  "*“kpt 

where  Lpb  and  Lpt  are  the  penetration  depths  of  the  bottom  and  top  mirrors  respectively 
which  have  been  derived  as 

tanh 2  (|k|-Lm  tb) 

2|k| 

where  LMt,bis  the  total  front  and  back  mirror  thicknesses,  and  k  is  the  coupling  constant 
which  is  given  by  approximately 

k  =  2  An/ X 

Here,  An  is  the  index  difference  between  the  two  layers  in  a  pair  of  the  mirror.  In  the 
standard  description  of  the  three-layer  waveguide,  the  TE  mode  is  normally  described  by 
a  function  within  the  guide  (asymmetric  or  symmetric)  and  by  evanescent  decay  into  the 
cladding  on  either  side  of  the  guide.  The  evanescent  decay  away  from  the  guide  is 
determined  by  the  refractive  index  in  the  region.  To  predict  the  behavior  of  the  vertical 
cavity  guide  we  assume  that  the  penetration  depth(s)  of  the  vertical  cavity  mirror(s) 
correspond  to  the  evanescent  decay  of  a  wave  if  it  were  propagating  in  the  waveguide 
formed  by  the  vertical  cavity  as  the  core  of  the  guide,  and  the  top  and  bottom  mirrors  as 
the  claddings.  Therefore,  by  using  Lpt,b  we  can  determine  the  effective  index  of  the  top 
and  bottom  mirror  regions  from  the  point  of  view  of  a  three-layer  waveguide.  The 
formula  given  for  Lpt,b  corresponds  to  a  vertically  traveling  wave  which  is  a  very  high 
order  mode  in  our  vertical  cavity  (VC)  waveguide,  i.e.  it  is  the  maximum  possible 
penetration.  However  in  the  fundamental  mode  the  V C  guide  will  propagate  at  some 
bounce  angle  and  the  penetration  will  be  less.  This  penetration  (and  thus  the  more 
accurate  indices)  will  be  determined  in  general  numerically  by  a  TMM  (transmission 
matrix  method)  solution  of  the  multi-layer  guide.  Using  these  indices  we  can  then  use  the 
model  developed  for  the  waveguide  diffraction  into  a  three  layer  waveguide  to  determine 
the  efficiency  of  the  second  order  diffraction  from  the  vertically  propagating  light  to  the 
guided  wave. 

For  the  Vertical  Cavity  device  we  have  designed,  the  mirror  pairs  on  the  top  are 
Si02  and  GaAs  corresponding  to  indices  of  1 .5  and  3.6.  The  bottom  mirror  is  comprised 
of  AlxOy  and  GaAs  corresponding  to  indices  of  1.6  and  3.6.  The  typical  design  will  be  7 
pairs  on  the  bottom  and  4-7  pairs  on  the  top.  Using  these  parameters  we  determine 

Lpt  =0.16  nm 

and  therefore  an  effective  index  above  the  core  of 
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Similar  calculations  for  the  lower  mirror  yield  3.05.  Then  the  diffraction  efficiency 
is  determined  from  our  grating  simulation  by  a  three-layer  waveguide  with  indices  of 
m= 2.96,  ri2  =3.48  and  m=  3.05.  The  index  of  the  core  region  (n2)  which  is  the  vertical 
cavity  itself,  is  determined  by  using  a  transmission  matrix  calculation  for  a  slab 
waveguide.  In  Fig.  5  we  plot  the  efficiency  of  the  diffraction  process  for  a  parallelogramic 
grating  in  the  VCSEL  waveguide  such  as  the  one  portrayed  in  Fig.l.  The  results  show  that 
we  can  expect  an  efficiency  of  0.0025  for  a  blaze  angle  of  38°  and  a  grating  etch  depth  of 
about  0.2  pm.  The  simulation  shows  an  efficiency  of  0.0033  at  an  angle  of  46°.  Such 
large  blaze  angles  may  be  difficult  to  achieve  with  ion  beam  etching  and  also,  for  a  depth 
of  0.2  pm  the  bottom  of  the  parallelogram  may  penetrate  back  to  undercut  the  top  of  the 
parallelogram  which  is  not  desirable.  We  thus  choose  an  angle  of  about  35°  which 
provides  a  value  of  diffraction  efficiency  of  about  0.00178.  For  this  calculation  we 
determine  that  the  ratio  of  the  forward  and  backward  diffracted  waves  is  600  or  28dB. 

This  is  an  important  consideration  for  rejection  of  reflected  waves  as  discussed  later. 


Tilt  angle  of  grating,  (degrees) 

Fig.  5  Diffraction  efficiency  versus  the  grating  angle  in  the  vertical  cavity  waveguide. 

It  is  noted  that  a  sharp  minimum  exists  in  the  efficiency  curve  at  a  certain  blaze  angle  of 
about  41°  .  This  effect  is  believed  to  result  when  the  angle  of  the  grating  coincides  with 
the  Brewster  angle  for  light  incident  normally  upon  the  waveguide.  It  is  currently  under 
investigation.  The  above  efficiencies  are  for  electric  fields.  Corresponding  power 
efficiencies  will  be  a  factor  of  2  larger 
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6.  DETAILED  DC  OPERATION  OF  THE  GCVCSEL 


From  the  above  description  one  can  see  that  the  light  generated  in  the  vertical 
cavity  at  z=0  will  be  diffracted  to  the  right  continuously.  The  diffracted  light  propagates 
in  the  form  of  a  guided  wave.  The  guided  wave  will  itself  be  diffracted  back  into  the 
cavity  at  each  position  z  and  will  be  amplified.  The  amplification  will  continue  until  the 
gain  is  saturated.  From  that  point  onwards,  the  saturated  vertical  cavity  output  will 
continue  to  add  to  the  waveguide  power  resulting  in  a  maximum  waveguide  output  power 
at  the  end  of  the  device  ,  z  =  L  . 

Consider  the  photon  density  at  any  position  z  in  the  cavity  of  thickness  dz  where 
the  laser  is  above  threshold.  We  have 

^cav  =  qL ”  -  Jth)  (1) 

where  Lx  is  the  quantum  well  width,  r\e  is  the  electrical  confinement  factor ,  J  and  JTH  are 
the  current  and  current  density  and  r'p  is  the  effective  photon  lifetime  for  the  vertical 


cavity  device  given  by 
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where  avc  is  the  vertical  cavity  loss  due  to  the  diffraction  grating.  There  will  also  be 
components  of  avc  which  are  due  to  free  carrier  absorption  and  parasitic  diffraction  but 
these  will  be  ignored  for  the  moment  to  focus  on  the  desired  effect.  In  the  final  formula 
we  may  simply  replace  ln(l/R)  by  Lcav  *apar  to  determine  its  effect  since  in  all  likelihood 
this  term  will  dominate  ln(l/R). 

The  power  output  from  the  laser  is  comprised  of  two  components  which  are 
written  as 


Pout  =  hvWdzvgt  fFcav 


(3) 


where  tf  is  the  power  transmissivity  of  the  output  port.  For  the  conventional  transmission 
through  the  mirror  (1)  is  substituted  into  (3)  with  tf  =ln(l/R)  to  give  the  power  increment 


dPmir  = 


hv 


Ve 


|  ^  Pcav^Vi 


-{j-JTH)Wdz 


(4) 


ln(l  /  R) 

There  is  another  power  component  which  is  the  power  diffracted  into  the  guide.  The 
transmissivity  for  this  mechanism  is  the  diffraction  efficiency  tf^rjdiff  so  we  have 


dP 


wg 


hvT__nsV, - (J-Jm)wdz 

qL-~ —ln(UR)  +  a,c 


(5) 


To  determine  avc  we  note  that  the  primary  optical  loss  in  the  vertical  cavity  by  design  will 
be  the  second  order  diffraction  by  the  grating  of  the  vertically  propagating  wave.  In  the 
development  of  diffraction  into  the  guide  in  this  program  [10]  from  a  normally 
propagating  wave  we  have  determined  the  diffraction  efficiency  rjdiff  as 
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where  rii  and  r|2  are  the  impedances  of  the  guide  and  the  incident  medium  (which  in  this 
case  are  the  same),  Ay  is  the  field  intensity  at  the  edge  of  the  guide  of  the  z  propagating 
wave  and  C0  is  the  field  intensity  of  the  incident  wave.  Also,  ag  is  the  growth  of  the 
guided  wave  due  to  diffraction  and  defy  is  the  effective  width  of  the  guide  which  includes 
evanescent  cladding  penetration.  Thus  we  can  say  that  the  power  diffracted  into  the  guide 
from  the  incident  wave  is  TidiffC0C0*  .  However  we  can  also  regard  the  grating  as  a  loss  in 
the  x  direction  for  the  wave  C0  and  we  can  define  a  loss  parameter  ax  by  the  statement 

c„c;(l-e--')  =  ^c„c;  (7) 


so  that  for  a  grating  thickness  of  x=t  we  have  (for  ax  x«l 

^  =r1d,ff 

Since  this  loss  occurs  only  over  the  grating  thickness  then  we  multiply  by  the 
confinement  factor  T  =  t  /  Lcav  to  obtain  an  effective  value  for  the  total  cavity  of 

avc  =  Vdiff  /  Lcav 

for  use  in  (5)  and  (2).  Therefore  the  waveguide  power  component  (5)  becomes 

dP  =— - y - 

wg  q  ln(l  /  R) V  ™ 

1  +  — - 

Vdiff 


(8) 

(9) 

(10) 


Since  the  concept  of  the  GC  VCSEL  is  to  maximize  the  waveguide  power  then  the  design 
of  the  laser  should  be  such  that  avc  is  the  dominant  loss  in  the  cavity  and  that  negligible 
power  escapes  through  the  two  mirrors  which  may  be  stated ,  by  using  (9)  in  (4),  as 

%#»ln(l/tf)  (11) 


From  (5)  it  is  seen  that  the  incremental  power  output  into  the  guide  is  just  the  total 
output  of  the  vertical  cavity  device  with  the  efficiency  of  the  diffraction  as  opposed  to 
mirror  transmission  as  in  (4).  According  to  (10)  the  power  in  the  guide  will  be 
proportional  to  the  guide  length.  As  the  wave  propagates  in  the  guide,  there  is  also 
diffraction  of  power  into  the  guide.  Basically  this  power  diffracts  from  the  wave  into  the 
cavity  and  then  back  into  the  guide  with  the  same  efficiency  as  the  cavity  flux.  However, 
this  photon  flux  can  amplify  the  cavity  flux  through  an  adjustment  of  Fcav  and  therefore  of 
(10).  In  general,  the  amplification  of  waveguide  power  depends  on  the  relative  size  of  the 
mirror  and  diffractive  loss.  If  the  diffractive  loss  is  very  large,  then  the  light  realizes 
single  pass  gain,  i.e.  it  passes  through  the  cavity  essentially  once  before  it  is  diffracted 
back  into  the  wave.  This  is  equivalent  to  the  linear  optical  amplifier  with  very  low 
reflectivity  in  which  case  the  majority  of  the  light  passes  through  the  amplifier  only  once 
as  a  traveling  wave.  In  these  cases  the  amplifier  always  works  below  the  threshold  of  the 
laser.  If  the  diffractive  losses  are  not  too  high,  then  the  light  may  make  many  oscillations 
in  the  cavity  before  diffraction  back  into  the  guide  and  this  situation  is  referred  to  as  the 
multipass  gain  case.  In  this  case,  the  laser  is  above  threshold  along  most ,  if  not  all  of  the 
length  of  the  laser.  To  obtain  maximum  laser  output,  this  is  the  situation  we  are  interested 
in  here.  For  this  kind  of  operation  (as  we  describe  elsewhere)  the  threshold  current  of  the 
laser  is  modified  .  We  can  therefore  write  an  equation  for  the  guide  power  which  is 
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where  in  this  case  (6)  is  modified  to 


F  = 

cav 


aFn 

jJzr3-) 


and  the  threshold  current  has  now  become  a  function  of  the  flux  injected  into  the  cavity. 
Here  Fg  is  the  waveguide  photon  density  given  by  Fg  =  Pwg/Avghv.  In  (12)  the  second 
term  is  essentially  (10)  and  the  first  term  is  the  net  diffraction  of  power  into  the  cavity.  It 
is  comprised  of  two  parts  which  are  first  the  diffraction  of  power  into  the  cavity  with 
efficiency  ag  per  unit  length  and  second  the  diffraction  of  those  same  photons  back  into 
the  guide  with  the  efficiency  we  derived  in  arriving  at  (10).  The  determination  of  the 
function  Jth^aFg/Fcav)  requires  a  modification  of  the  photon  rate  equation  from 


~  cav  r  cavT  =  L  cav 


aFgvg 


The  equation  on  the  right  is  the  conventional  photon  loss  equation  where  the  right  hand 
side  represents  photon  loss  from  the  system.  The  equation  on  the  right  is  the  photon  loss 
equation  when  an  optical  input  is  fed  to  the  device.  It  is  clear  that  the  input  photon  term 
reduces  the  loss  and  therefore  reduces  the  condition  for  threshold  and  therefore  the 
threshold  current.  Specifically ,  the  K  parameter  [9]  which  is  used  to  determine  threshold, 
contains  x'p  and  from  (14) ,  the  photon  lifetime  is  modified  to 


=  Cd 


«gFgvg*P, 


i.e.,  the  effective  photon  lifetime  in  the  cavity  increases  with  the  input  signal  but 
decreases  with  the  cavity  flux  itself. 

The  power  output  of  the  laser  is  determined  by  integration  of  (12)  and  is  shown  as 
a  function  of  laser  length  L  in  Fig.6.  If  the  effect  of  optical  gain  in  the  cavity  were  absent 
so  that  Jth  in  (12)  is  a  constant  then  the  solution  to  (12)  is 


p-<L> - p  ■  j- j"  °6) 

which  is  equivalent  to  the  output  of  a  vertical  cavity  device  with  a  width  of  W  ,  a  length 
of  L=  l/ag  and  an  efficiency  of  r|=  pdiff  Pe  /{ln(l/R)+Lcav*apai).  Equivalently,  we  could 
say  the  efficiency  was  r|eand  the  effective  length  was  t|diff  /ag  { ln(  1  /R)+Lcav* otpar) .  In  the 
figure,  the  parameters  are  taken  from  a  typical  grating  efficiency  analysis  and  the 
parameters  of  ridiff  =0.001  and  ag=50  are  used  with  a  cavity  designed  for  R=0.9999. 

Then  the  effective  length  is  0.7  cm  .  This  may  be  stated  alternatively  by  saying  the  device 
has  an  efficiency  of  63%  for  a  length  of  L  =  0.7  cm.  Desirable  parameters  are  therefore 
r|=0.001  and  ag=50  which  can  be  achieved  with  the  dielectric  combination  SiCVGaAs 
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and  a  grating  etch  depth  (thickness)  of  1000  angstrom.  From  our  calculations  for  an 
asymmetrical  grating  the  optimum  blaze  angle  is  about  30°.  With  these  numbers  ,  the 
laser  output  power  for  J-Jth  =  834A/cm2  and  a  waveguide  width  of  10pm  is  about  0.24W. 


Device  length,  (cm) 

Fig.  6  Waveguide  photon  density  versus  device  length 


Also  shown  on  the  plot  is  the  output  power  versus  z  when  the  dependence  of  Jth 
on  input  power  is  included  in  (12).  It  is  evident  that  the  gain  internal  to  the  device 
increases  the  output  power  by  a  factor  of.  The  overall  efficiency  of  the  laser  is  a  function 
of  z  and  therefore  of  the  actual  device  length  and  from  (16)  we  obtain 


ri  = 


Wle 

agL/«(l/R) 


|l-exp 


-a. 


/«(  i/r) 

ri  dS+ln(l/R) 


One  interesting  feature  is  the  tradeoff  between  length  and  overall  efficiency.  For  small 
z=L,  the  efficiency  is  close  to  100%.  As  the  length  is  increased  the  efficiency  drops  so 
that  at  L«0.7cmthe  efficiency  has  dropped  to  about  45%.  Note  however  that  we  cannot 
use  the  prefactor  of  J-Jth  in  (16)  as  the  efficiency  since  this  result  implies  that  the  power  is 
independent  of  z  (i.e.  as  the  length  is  reduced  to  zero,  the  exponential  term  reduces  the 
output  power  to  zero).  This  is  clearly  a  design  tradeoff  in  which  the  total  power  and 
efficiency  must  be  traded  off  by  the  appropriate  choice  of  length. 

The  operation  of  the  device  as  an  optical  amplifier  will  not  be  considered  in  detail 
here  because  the  laser  is  our  main  objective.  However,  it  is  interesting  to  note  that  as  long 
as  the  rjdiff  is  large  enough  to  prevent  lasing  it  may  still  be  quite  small,  i.e.  it  may  easily 
be  around  10'3.  The  gain  is  obtained  because  for  every  increment  of  propagation  distance 
in  the  z  direction  the  wave  traverses  the  vertical  cavity  Q  times  where  Q  is  the  quality 
factor.  Since  Q=  (Vq/v f)3  where  3  is  the  finesse  given  by  approximately  3=7t/at=  nh idiff 
and  v0  and  vp  are  the  resonator  frequency  and  free  spectral  range  respectively  then  the 
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light  may  traverse  a  distance  of  QLcav  vertically  in  a  propagation  distance  of  dz  in  the 
guide.  Estimating  that  1-1/e  of  the  wave  is  diffracted  into  the  cavity  in  the  distance  l/ag 
then  the  effective  distance  traveled  by  the  wave  is  about  (v0/vf)  nlr\^.L  cav  •CtgLguide 
during  which  it  receives  optical  gain.  Using  typical  numbers  this  is  a  distance  of  15cm.. 
This  indicates  the  feasibility  of  implementing  an  optical  amplifier  on  a  chip  with 
monolithic  integration  which  has  all  the  benefits  of  the  fiber  amplifier.  This  points  to  the 
direction  of  implementing  an  integrated  transmission  system. 

There  are  other  benefits  to  be  realized  from  this  grating  coupled  arrangement 
which  are: 

1)  the  polarization  of  the  vertical  cavity  wave  is  forced  to  coincide  with  that  of 
the  waveguide.  The  locking  of  the  polarization  occurs  because  the  waveguide  mode  is 
injected  into  the  cavity  at  each  point  and  creates  some  level  of  stimulated  emission.  The 
emission  in  the  cavity  reproduces  the  same  polarization.  The  random  nature  of  the 
polarization  of  the  vertical  cavity  laser  output  has  now  been  eliminated.  The  output  of  the 
waveguide  will  be  TEo  normally  because  it  is  most  easily  excited  in  the  guide.  Earlier 
simulations  of  Lee  and  Streifer  [11]  have  shown  in  general  that  the  diffraction  efficiency 
of  the  TMo  mode  is  at  least  10  times  smaller  than  the  TE0  mode,  which  thus  becomes  the 
principle  supported  mode.  Because  of  the  diffraction  of  the  light  from  the  guide  into  the 
cavity  there  will  be  a  stabilizing  effect  and  an  extended  range  of  single  mode  output 
power.  For  a  conventional  VCSEL,  the  output  is  single  mode  until  a  certain  level  of 
power  is  achieved  and  then  multimode  behavior  is  observed.  In  the  grating  coupled  laser 
the  onset  of  higher  order  VCSEL  modes  will  require  a  much  higher  level  of  laser  power 
to  occur  because  of  the  stabilizing  effect  of  the  waveguide  injected  energy.  In  fact,  the 
GCVCSEL  will  remain  single  mode  until  a  higher  order  mode  in  the  waveguide  matches 
to  a  higher  order  mode  in  the  VCSEL. 

2)  The  mode  supression  ratio  will  be  characteristic  of  the  vertical  cavity  laser. 
Because  of  the  single  mode  nature  of  the  VCSEL,  the  tendency  for  mode  partition  noise 
will  be  much  reduced.  The  RIN  noise  should  be  typical  of  the  VCSEL 

3)  One  of  the  problems  of  edge  emitting  and  vertical  cavity  lasers  is  the  effect  of 
back  reflections  into  the  laser  when  coupling  to  a  fiber.  In  the  grating  coupled  VCSEL, 
this  problem  is  addressed  in  a  unique  way  as  shown  in  Fig.2.  Light  is  coupled  from  the 
laser  to  the  waveguide  and  is  transported  to  the  chip  edge.  Reflections  travel  back  to  the 
laser  and  a  standing  wave  will  be  set  up  which  may  be  at  a  different  frequency.  Normally 
this  energy  will  re-enter  the  cavity  and  may  destabilize  the  laser.  In  the  grating  coupled 
laser,  which  we  have  modeled  as  an  equivalent  three-level  waveguide,  the  asymmetry  of 
the  grating  causes  a  dominant  fraction  of  the  power  to  be  diffracted  into  the  guide  for 
light  traveling  towards  the  chip  edge  (direction  of  laser  emission)  and  it  also  causes  a 
dominant  fraction  of  the  light  to  be  diffracted  out  of  the  waveguide  for  light  approaching 
the  laser  from  the  edge  of  the  chip.  In  essence  then,  the  grating  asymmetry  acts  to  reject 
the  light  re-entering  the  laser  with  a  ratio  of  diffraction  up  to  diffraction  down  which  can 
be  up  to  28dB  according  to  the  simulation.  The  grating  thus  acts  as  an  optical  isolator. 

4)  The  GCVCSEL  output  is  into  a  waveguide  on  the  chip.  The  waveguide 
transports  the  light  to  the  chip  edge.  Output  coupling  of  the  laser  light  to  a  fiber  occurs 
from  the  waveguide  to  the  fiber.  One  of  the  major  problems  in  manufacturing  involves 
attaching  fiber  to  the  chip  and  much  effort  has  gone  into  the  development  of  spot  size 
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transformers  which  are  ways  to  match  the  waveguide  mode  size  to  the  fiber  aperture.  The 
waveguide  output  is  easily  tapered  in  the  lateral  dimension  to  achieve  this  matching.  In 
the  vertical  dimension  the  mode  size  can  be  enlarged  considerably  by  the  use  of  the  ion 
implant  to  disorder  the  waveguide.  As  we  point  out  in  the  process  flow  sequence  (section 
8)  the  waveguide  core  is  implanted  with  Si  during  the  formation  of  the  source-drain 
regions.  This  implant  after  RTA  serves  to  disorder  the  bandgap  of  the  quantum  well 
layers  which  reduces  the  index.  The  net  effect  is  to  move  the  core  index  closer  to  the 
effective  index  of  the  upper  and  lower  quantum  well  regions  which  causes  the  expansion 
of  the  mode  into  these  regions.  Judicious  use  of  this  implant  and  anneal  step  can  be  used 
to  optimize  the  mode  shape  to  increase  the  tolerances  for  the  mode  to  fiber  coupling 
problem. 

Another  way  to  widen  the  mode  shape  vertically  is  to  remove  layers  of  cladding 
from  the  waveguide  as  the  chip  edge  is  approached.  In  the  device  area  there  are  many 
dielectric  pairs  to  eliminate  emission  through  the  upper  mirror.  This  maintains  tight  mode 
confinement.  Along  the  waveguide  these  layers  can  be  selectively  removed  to  widen  the 
mode  shape. 
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7.  HIGH  SPEED  OPERATION  OF  THE  GCVSEL 


A  primary  motivation  for  this  design  of  laser  is  to  achieve  very  high  speed 
performance.  Within  the  inversion  channel  technology  the  HFET  laser  is  ideally  suited  to 
achieve  the  necessary  band  widths.  The  HFET  laser  is  a  laterally  injected  laser  in  which 
the  gate/emitter  and  source  injecting  terminals  lay  next  to  each  other  along  the  surface 
and  are  therefore  configured  as  adjacent  electrodes  as  illustrated  in  Fig.2.  Also,  in  the 
vertical  cavity  arrangement  as  illustrated  in  Fig.  lor  2,  the  distance  from  the  quantum 
wells  to  the  semiconductor  surface  can  be  54  wavelength  as  in  most  VCSEL  designs  and 
since  the  source  drain  implant  is  performed  from  approximately  1 000  A  above  the 
quantum  well  to  achieve  minimum  channel  access  resistance,  then  the  positive  and 
negative  electrodes  form  an  almost  perfect  coplanar  transmission  line.  The  vertical  cavity 
of  the  HFET  laser  structure  is  formed  by  the  deposition  of  the  dielectric  stack  over  the 
electrode  structure.  These  dielectric  layers  form  the  top  dielectric  cover  for  the  coplanar 
line  and  are  therefore  important  in  determining  Z0.  When  implemented  in  this  form  the 
HFET  laser  implements  almost  perfectly  a  traveling  wave  laser.  This  traveling  wave 
property  extends  the  bandwidth  of  the  laser  to  the  limit  imposed  by  the  internal  parasitic 
resistance  and  time  constant  because  the  coplanar  electrical  transmission  line  can  be 
terminated  on  the  chip  by  a  transistor  adjusted  in  impedance  to  match  the  Z0  of  the  line. 
Then  if  the  transmission  line  phase  velocity  and  the  waveguide  group  velocity  are 
matched ,  one  could  introduce  a  signal  on  the  coplanar  line  which  would  propagate  at  the 
same  speed  as  the  photon  pulse  which  was  being  produced  by  the  continuous  injection  of 
charge  into  the  channel  from  the  pulse  propagating  on  the  line.  Suppose  there  was  no 
dispersion  of  the  electrical  pulse  and  no  dispersion  of  the  optical  pulse.  If  there  was 
negligible  delay  in  producing  the  photon  (i.e.  negligible  transit  time  into  the  channel)  and 
negligible  delay  in  converting  the  electron  to  a  photon  (very  fast  laser  internal  response) 
then  the  laser  would  have  an  infinite  bandwidth  since  then  a  delta  function  of  charge 
introduced  to  the  line  would  produce  (after  the  transit  time  of  the  waveguide)  a  delta 
function  of  photons.  The  extent  to  which  these  conditions  are  not  met  determines  the 
maximum  laser  speed  and  generally  it  is  the  internal  parasitic  RC  constant  which 
becomes  the  limiting  factor(as  we  show  in  the  detector  section  7  this  parasitic  bandwidth 
is  about  100GHz). 

To  determine  further  the  dynamic  operation  consider  a  pulse  introduced  on  to  the 
coplanar  line  at  x=0.  Because  of  the  almost  identical  structure  of  the  optical  waveguide 
and  the  transmission  line  in  terms  of  the  materials  both  above  and  below  the 
semiconductor  surface,  the  optical  group  velocity  and  the  electrical  phase  velocity  will  be 
almost  matched.  As  the  pulse  propagates  on  the  line,  it  will  continuously  decrease  in 
value  as  charge  is  injected  laterally  at  each  position  x  into  the  laser  active  section.  Now 
the  effect  of  the  injected  charge  will  be  a  maximum  if  the  charge  is  injected  at  precisely 
the  point  of  the  maximum  photon  density.  Also  the  charge  impulse  initially  begins  to 
create  a  photon  impulse  starting  at  x=0.  Then  if  the  photon  pulse  and  the  charge  pulse 
move  at  the  same  velocity  (i.e.  the  matched  condition),  it  follows  that  as  the  charge 
impulse  travels  on  the  transmission  line  it  will  always  be  injecting  its  incremental  charge 
into  the  laser  at  the  position  which  corresponds  to  the  maximum  photon  density  in  the 
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waveguide.  If  this  is  the  case,  the  photon  pulse  will  grow  the  most  efficiently  because  the 
number  of  photons  produced  by  the  injecting  charge  is  proportional  to  the  injected  charge 
and  the  photon  density,  i.e.  the  product  of  the  electron  and  photon  densities.  This 
behavior  is  illustrated  schematically  in  Fig.7  which  shows  the  case  at  some  position  x 
along  the  waveguide  the  situations  when  the  electrical  pulse  travels  slower,  faster  and  at 
exactly  the  same  speed  as  the  optical  pulse. 
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Fig.  7  Effective  velocity  matching  on  the  optical  output  response 

Ideally  when  they  travel  at  the  same  speed,  the  photon  pulse  is  a  maximum  height 
because  the  injected  charge  x  photon  density  product  has  been  a  maximum  at  all  points 
during  the  propagation.  Also,  the  pulse  width  will  be  the  narrowest  because  ,  assuming 
that  the  internal  laser  dynamics  are  sufficiently  fast ,  the  pulse  width  of  the  photons  will 
be  as  narrow  as  the  charge  pulse.  If  there  is  negligible  dispersion  on  the  transmission  line 
then  the  charge  pulse  and  thus  the  photon  pulse  will  retain  its  original  delta  function  form, 
i.e.  the  photon  pulse  will  remain  essentially  as  fast  as  the  input  electrical  pulse. 

A  complete  analysis  of  the  laser  dynamic  performance  is  exceedingly  complicated 
if  one  attempts  to  illustrate  the  effects  of  the  velocity  matching  or  lack  thereof.  This  is 
because  of  several  properties  unique  to  the  lasing  mechanism  which  are 

•the  photon  response  is  proportional  to  the  product  of  the  injected  carriers  and  the 
injected  photons.  Therefore  in  the  laser  cavity,  at  any  position  along  the  guide  we  need  to 
consider  a  mixture  of  the  dispersed  pulse  traveling  on  the  transmission  line  and  the 
dispersed  pulse  traveling  in  the  waveguide  and  to  represent  their  individual  time 
dependencies. 

•in  the  laser  there  is  an  amplifier  effect  which  means  photons  produce  additional 
photons  without  the  need  for  additional  electrons.  This  may  be  considered  an  optical  to 
optical  effect  as  opposed  to  the  electrical  to  optical  effect  of  the  carriers  injected  from  the 
transmission  line.  Therefore  we  need  to  consider  the  dynamics  of  the  optical  amplifier 
effect  in  the  waveguide  itself  so  that  the  photon  pulse  is  growing  not  only  due  to  the 
injection  from  the  line  but  also  due  to  the  internal  optical  gain.  And  as  we  mentioned 
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earlier,  this  optical  gain  can  be  in  one  of  two  forms  depending  upon  how  far  above  or 
below  threshold  the  laser  is  biased.  This  situation  is  to  be  contrasted  with  the  detector  in 
which  electrons  are  produced  (i.e.  the  reverse  process)  only  by  the  conversion  of  the  pulse 
photons. 

•proper  analysis  requires  use  of  the  actual  differential  equations  representing  the 
electron  and  photon  populations  but  represented  in  the  form  of  traveling  waves.  These 
equations  have  not  yet  been  solved  self-consistently  in  the  literature  for  the  conventional 
diode  laser  without  any  traveling  wave  effects. 

Therefore  more  work  is  needed  to  obtain  a  clear  analysis.  However ,  the  physical 
picture  described  above  is  clear  which  is  that  if  the  waveguide  and  transmission 
dispersions  can  be  minimized  and  the  internal  parasitic  RC  constant  can  be  reduced  to  5- 
1  Ops,  then  velocity  matching  of  the  HFET  transmission  line  and  the  vertical  cavity 
waveguide  can  enable  speeds  in  excess  of  100GHz. 

The  analysis  of  the  traveling  wave  detector  is  much  easier  and  reinforces  all  of  the 
conclusions  made  here  for  the  laser.  It  is  presented  in  the  next  section. 
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8.  THE  HFET  VERTICAL  CAVITY  TRAVELING  WAVE  DETECTOR 


8.1  General  Design 

To  meet  the  transmission  requirements  for  detection  requires  innovation  also  in 
the  detector.  Traditional  detector  geometries  available  to  the  designer  include  the 
classical  vertically  illuminated  PIN  structure  or  the  waveguide  structure  which  is  edge 
illuminated .  The  inherent  bandwidth-efficiency  trade-off  in  the  vertical  device  between 
absorption  efficiency  in  the  active  layer  and  carrier  transit  time  through  it  have  been 
addressed  significantly  by  the  implementation  of  resonant  cavity  enhanced  (RCE) 
structures  wherein  multiple  optical  passes  are  used  to  obtain  almost  total  absorption  in  a 
very  narrow  active  layer ,  i.e.  in  a  quantum  well.  The  drawback  of  the  RCE  is  the  narrow 
optical  bandwidth  which  results  from  the  high  finesse  of  the  cavity.  Nevertheless,  for 
many  of  the  emerging  practical  system  architectures  in  which  wavelength  division 
multiplexing  is  employed,  the  wavelength  selectivity  offered  by  the  cavity  is  exactly  what 
is  required  to  perform  optical  to  electronic  demultiplexing.  The  difficulty  of  course  is  to 
devise  a  means  to  inject  the  light  into  the  cavity  since  the  high  reflectivity  of  the  cavity 
will  reject  the  optical  input  unless  it  is  exactly  positioned  within  the  narrow  bandwidth 
corresponding  to  the  position  of  the  cavity  optical  mode. 

In  the  waveguide  configuration,  it  is  typical  to  use  a  double  heterostructure 
semiconductor  structure  with  a  quantum  well  active  region.  The  thickness  of  the  active 
layer  can  also  be  minimal  since  the  absorption  occurs  along  the  length  of  the  guide  which 
can  be  designed  for  complete  absorption.  The  transit  time  limitation  and  the  RC 
limitation  can  be  essentially  the  same  as  that  of  the  vertically  illuminated  RCE  device. 
The  efficiency  of  the  device  does  suffer,  however,  from  the  problem  of  poor  input 
coupling  since  the  numerical  aperture  of  the  typical  semiconductor  waveguide  is  small 
and  the  mode  is  not  well  matched  to  the  typical  mode  of  the  optical  fiber.  These  devices 
have  been  reported  with  50  GHZ  bandwidth  and  40%  external  quantum  efficiency  and 
also  68%  external  quantum  efficiency  with  >50GHZ  bandwidth. 

Both  of  the  above  approaches  behave  as  lumped  element  devices  electrically  and 
consequently,  the  reported  bandwidths  are  approximately  limited  by  the  load  resistance 
and  the  device  capacitance.  To  alleviate  the  lumped  element  limitation  and  thereby 
improve  bandwidth,  traveling  wave  structures  have  been  proposed.  These  devices 
basically  combine  the  optical  waveguide  configuration  with  an  electrical  transmission 
line.  This  approach,  which  was  introduced  initially  to  enhance  the  speed  of  lithium 
niobate  waveguide  modulators  has  the  potential  to  substantially  reduce  the  effect  of  the 
RC  time  constants  imposed  by  external  circuits.  The  advantage  is  obtained  because  the 
load  impedance  presented  to  the  detector  does  not  combine  with  the  internal  impedance 
of  the  detector  to  establish  a  time  constant  if  the  load  is  matched  to  the  characteristic 
impedance  of  the  line.  The  temporal  response  is  further  improved  if  the  transmission  line 
energy  velocity  and  the  waveguide  optical  group  velocity  are  matched,  i.e.  the  condition 
of  zero  velocity  dispersion. 

In  practice,  impedance  matching  and  velocity  matching  are  very  difficult. 
Furthermore,  the  TWPD  still  suffers  from  the  problem  of  poor  input  coupling  since  the 
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basic  semiconductor  waveguide  structure  has  remained  unchanged.  In  fact  the  waveguide 
must  be  lpm  or  less  in  width  to  realize  reasonable  transmission  line  parameters  and  this 
exacerbates  the  coupling  problem. 

Another  limitation  with  the  existing  traveling  wave,  waveguide  or  RCE  designs  is 
the  problem  of  integration  with  electronic  devices.  A  traveling  wave  device  intended  for 
1 00  GHz  operation  must  feed  photocurrent  to  an  FET  or  bipolar  front  end  in  a  matched 
configuration  to  avoid  reflections.  It  is  not  practical  or  cost  effective  to  do  this  with 
hybrid 

connections  except  for  specialized  applications  and  it  thus  becomes  essential  to  have  an 
integrated  approach.  Ideally,  the  integration  should  encompass  the  detector,  the  laser,  and 
the  electronic  amplifier. 

Since  we  intend  to  use  the  Inversion  channel  technology  approach  we  then 
propose  the  HFET  resonant  cavity  enhanced  detector  as  a  novel  means  to  address  these 
problems  since  it  combines  the  virtues  of  the  vertical  cavity  RCE  structure  and  the  TW 
concept  into  a  single  device.  This  approach  has  several  advantages  which  include  1) 
optimized  mode  matching  to  an  optical  fiber  to  achieve  improved  input  coupling,  2) 
optimized  impedance  matching  to  reduce  reflections,  3)  optimized  velocity  matching  to 
reduce  pulse  dispersion,  and  4)  integrated  circuit  compatibility  to  achieve  a  low  cost  high 
performance  package. 

In  the  ICT  the  laser  and  detector  are  physically  the  same  structure.  Therefore  the 
Figs.  1  and  2  also  apply  here.  These  figures  show  the  traveling  wave  version  of  the 
HFET  laser  or  detector  since  physically  these  are  the  same  structure.  To  differentiate 
between  the  laser  and  the  detector  operation  requires  the  device  cross-section  in  Fig.8. 
which  shows  the  appropriate  current  flows. 

9 


Laser 


Detector 

Fig.  8  Current  flow  diagram  for  laser  and  detector 
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This  figure  indicates  that  the  collection  of  carriers  in  the  devices  is  lateral  for  the 
electrons  and  vertical  for  the  holes.  If  the  collector  is  contacted  with  a  negative  bias  then 
the  holes  are  removed  downward  to  the  collector.  On  the  other  hand  if  the  gate  is 
contacted  with  a  negative  bias,  the  hole  current  flows  upward  to  the  gate  electrode.  It  is 
noted  that  the  collector  operation  can  be  achieved  with  zero  collector  to  source  bias  since 
the  barrier  at  zero  bias  is  low  whereas  the  gate  operation  requires  a  bias  to  lower  the 
channel  to  gate  barrier  for  holes  preferentially  relative  to  the  collector.  The  electrons  flow 
within  the  inversion  channel  and  are  collected  at  the  source  nodes.  The  electrons  flow  to 
both  source  nodes  so  that  the  worst  case  transit  time  corresponds  to  1/2  of  the  channel 
length  and  the  average  effective  delay  corresponds  to  a  transit  distance  of  L/4.  Note  that 
the  vertical  cavity  geometry  in  Fig.2  will  typically  admit  light  through  an  opening  in  the 
gate  electrode  so  that  the  conduction  of  holes  to  the  gate  is  by  the  2D  funneling 
mechanism  which  is  typically  used  in  the  VCSEL  structure  to  provide  current  injection 
into  the  same  region  as  the  optical  emission.  No  conduction  is  required  through  the  mirror 
structure  however  since  the  mirror  is  applied  in  the  form  of  a  dielectric  stack  at  the  end  of 
the  fabrication  sequence..  For  the  case  of  the  waveguide  configuration  in  Fig.2,  the  gate 
electrode  is  continuous  since  the  wave  is  guided  under  the  top  contact  and  sufficiently 
removed  to  avoid  excessive  loss. 

The  description  of  the  technology  changes  to  accommodate  the  detector  are 
essentially  the  same  as  those  of  the  laser  above.  The  metal  contact  is  opened,  as  in  the 
conventional  vertical  cavity  device  ,  to  allow  the  passage  of  the  light  in  the  vertical 
direction.  However  here  the  light  is  propagating  both  horizontally  and  vertically.  The 
dielectric  mirror  of  the  vertical  cavity  device  is  used  here  in  a  multi-functional  role.  First, 
the  semiconductor  growth  between  the  quantum  well  active  region  and  the  dielectric 
stack  is  very  thin,  i.e.  approximately  A./2n ,  where  n  is  the  average  material  index,  in  order 
to  move  the  gate  contact  as  close  as  possible  to  the  inversion  channel.  This  is  necessary  to 
form  a  high  frequency  transistor  structure  .  However  the  optical  mode  is  still  efficiently 
guided  because  the  dielectric  stack  provides  efficient  guiding  of  the  mode.  The  mode  is 
centered  in  the  quantum  wells  but  extends  well  above  the  gate  metal  contact.  Second,  the 
dielectric  stack  also  forms  a  high  finesse  cavity  in  the  vertical  direction.  Therefore  the 
number  of  pairs  in  the  stack  and  the  index  difference  are  selected  to  achieve  the  desired 
finesse  for  the  resonant  cavity.  Third,  the  multiple  layer  structure  of  the  cavity  mirror 
allows  a  large  optical  mode  to  be  supported  for  propagation  in  the  horizontal  direction. 
The  enlarged  mode  size  increases  the  mode  near-field  pattern  and  decreases  the  far-field 
angles  which  facilitates  easier  coupling  to  a  fiber.  The  input  light  is  edge  coupled  to  a 
waveguide  at  the  chip  edge  with  a  cross-section  designed  for  optimum  coupling.  The 
waveguide  guides  the  light  to  the  detector  where  it  propagates  along  the  device 
waveguide  as  shown  in  Fig.l. 

As  for  the  laser  the  key  element  to  enable  the  combined  waveguide,  vertical  cavity 
operation  is  a  grating  which  is  created  in  the  first  layer  of  the  dielectric  stack.  This  is  a 
second  order  grating  which  diffracts  a  portion  of  the  input  wave  into  the  vertical  cavity  at 
each  position  along  the  guide.  Since  the  light  propagating  vertically  in  the  cavity  is 
absorbed  essentially  completely  due  to  the  resonant  enhancement  effect  then  the 
limitation  on  the  detector  efficiency,  i.e.  the  length  of  detector  required,  is  determined  by 
the  grating  parameter  ag  .  Because  the  grating  index  change  can  be  made  fairly  large  by 
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suitable  design  of  the  dielectric  stack  then  ag  can  be  large  and  the  waveguide  length  can 
be  optimized  for  high  speed  operation. 

One  of  the  obvious  features  of  this  design  is  that  since  the  light  is  transferred  to 
the  resonant  enhanced  structure  via  the  grating  which  passes  below  the  dielectric  stack 
then  the  problem  of  injecting  light  into  the  cavity  through  a  dielectric  mirror  with  a  very 
narrow  modal  bandwidth  has  been  eliminated.  Clearly,  the  light  will  not  be  absorbed  if  its 
frequency  does  not  align  with  that  of  the  mode.  However,  this  is  a  desirable  property, 
since  it  allows  the  cavity  to  select  any  particular  frequency  as  would  be  required  in  a 
WDM  application,  i.e.  the  light  is  always  coupled  into  the  cavity  with  the  cavity  mode 
providing  the  filtering  or  demultiplexing  operation.  Those  frequencies  which  are  not 
selected  are  simply  coupled  back  into  the  waveguide  by  the  grating  and  are  propagated 
out  of  the  device. 

The  HFET  TW  detector  illustrated  in  Fig.2  uses  the  gate  contact  to  collect  the 
hole  current  in  order  that  it  may  form  the  signal  line  for  the  transmission  line.  The  other 
side  of  the  line,  i.e.  the  ground  line,  is  formed  by  the  metal  electrodes  on  the  source 
junctions.  Due  to  the  aspect  ratio  as  indicated  in  Fig.3,  which  shows  the  two  metal 
electrodes  at  approximately  the  same  elevation,  (within  1500-2000  A),  the  source  and  gate 
electrodes  behave  as  a  coplanar  transmission  line.  Therefore,  as  photo-charge  is  injected 
into  the  source  and  gate  electrodes,  these  perturbations  form  a  traveling  wave  on  the 
coplanar  line  which  propagates  to  the  far  end  of  the  device(x  =  L  in  Fig.3).If  the  group 
velocity  of  the  transmission  line  wave  and  the  group  velocity  of  the  optical  wave  are 
matched  then  the  photo-charge  is  added  synchronously  to  the  electrical  wave  and  there 
will  be  no  distortion  of  the  pulse  due  to  velocity  mismatch.  If  the  transit  time  of  holes  to 
the  gate  from  the  channel  and  of  electrons  to  the  source  contact  from  the  channel  are 
equal  and  if  there  is  no  dispersion  on  the  transmission  line  or  the  waveguide  then  the 
input  pulse  shape  will  be  preserved.  Finally  if  the  transmission  line  is  exactly  matched  in 
its  characteristic  impedance,  there  will  be  no  reflection  and  the  pulse  will  be  perfectly 
replicated,  i.e.  the  impulse  response  of  the  detector  would  be  infinite.  These  conditions 
can  never  be  realized  perfectly  and  it  is  the  deviation  of  these  criteria  from  the  ideal 
which  determines  the  actual  detector  bandwidth. 

The  virtue  of  the  traveling  wave  concept  is  that  the  impedance  of  the  circuit 
following  the  detector  does  not  combine  with  the  detector  impedance  to  produce  delay  if 
ideal  line  matching  is  achieved.  However  parasitic  RC  delays  which  are  intrinsic  to  the 
device  are  still  an  issue.  These  are  (Fig.2)  the  time  constant  associated  with  transferring 
the  electron  and  hole  from  the  absorption  region  to  the  source  metal  and  gate  metal 
respectively.  To  minimize  these  parasitic  effects,  the  transmission  line  metal  to  metal 
spacing  should  be  reduced  and  the  channel  doping  increased  as  much  as  possible  by 
aggressive  technology  scaling.  The  other  fundamental  delay  is  the  transit  time  of  electrons 
in  the  absorbing  FET  channel.  For  a  channel  length  of  L  =1  pm,  this  time  corresponds  to 
an  average  of  L/4vs  For  a  vs  «2*107  cm/sec  at  low  fields,  this  corresponds  to  a  delay  of 
about  0.15  ps.  From  the  point  of  view  of  efficiency,  the  attenuation  of  the  transmission 
line  ae  must  now  be  added  to  the  losses  in  the  resonant  cavity  enhanced  detection  process. 
This  attenuation  increases  with  frequency  with  dependencies  on  the  source  and  gate  metal 
skin  effect  depths  and  the  frequency  dependent  propagation  constant.  In  this  respect,  the 
bandwidth  limitation  of  the  lumped  element  configuration  has  been  exchanged  for  a 
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conversion  efficiency  limitation  resulting  from  the  attenuation  or  loss  of  the  traveling 
photo-induced  current. 

8.2  Dynamic  Response 

The  electrical  model  of  the  transmission  line  within  the  detector  is  shown  in  Fig.9. 
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Fig.  9  Circuit  model  of  transmission  line  for  traveling  wave  optoelectronic  devices 

The  electrical  losses  in  the  transmission  are  represented  by  R  and  G.  The  series  R 
of  the  transmission  line  is  determined  by  the  combined  skin  effect  resistance  of  the  source 
and  gate  conductors  which  may  be  expressed 


R  =  Rs  +  Rg  - 
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where  co  =  2nf  is  the  frequency,  ws  and  wg  are  the  source  and  gate  widths,  p  is  the 
permittivity  and  as  ,  crg  are  the  metal  conductivities.  The  parallel  conductance  represents 

essentially  dielectric  loss  mechanisms  in  addition  to  dark  current  leakage.  These 
parasitics  combine  to  give  a  total  loss  parameter  a  such  that  the  wave  on  the  transmission 
line  propagates  as 

E{x)  =  E0e'cae-Jfix  (18) 


where  both  a  and  p  are  frequency  dependent,  as  outlined  in  standard  textbooks.  The  loss 
may  be  written  as  the  sum  of  the  series  and  parallel  contributions  as 

ae=ac+ad  (19) 


and  for  the  unbalanced  coplanar  line,  these  are  expressed 
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and 
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where  £(*:,),  £(£[),£(£  ),£(£')  are  the  complete  elliptic  integrals  of  the  first  kind,  and 


k\  =  -J\-k2  .  The  effective  index  is  also  frequency  dependent  and  is  expressed  [12]  as 

=  '  •*  -  PD 

where 


1  +  f  /  fTE  ) 


A  =  exp 


0.54 -0.64  log! 

2  a 


2  a 
~h) 


r 


+  0.015 


log 


f2a^ 


-  0.86  log[  —  j  +  0.54| 


iog| 


KhJJ 

h  ) 


log 


2  a  ^ 


b-aJ 


+  0.43 


f TE 


^0 


The  impulse  response  of  the  TW  detector  has  been  stated  as 
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which  represents  both  forward  and  backward  traveling  waves,  where  E0  is  the  delta 
function  input  impulse  of  energy,  a  is  the  absorption  coefficient  of  the  active  layer,  T  is 
the  optical  confinement  factor,  v0  is  the  optical  velocity,  veis  the  electrical  energy  ( or 
group )  velocity,  /  is  the  detector  length,  and  t  is  the  time  as  measured  from  the  time  of 
incidence  of  energy  to  the  detector,  77,  is  the  internal  quantum  efficiency,  q  is  the 
electronic  charge,  hv  is  the  photon  energy.  In  this  section  we  will  extend  this  result  to 
demonstrate  the  importance  of  the  electrical  loss  on  both  the  detector  efficiency  and 
velocity  mismatch  which  then  establishes  the  bandwidth  limit.  We  then  utilize  the  result 
to  determine  the  frequency  dependent  response  to  a  realistic  pulse  input. 


8.2.1  Impulse  Response  and  Small  Signal  Transfer  Function 

The  optical  pulse  enters  the  waveguide  as  shown  in  Fig.  10  at  x  =  0  and  continuously 
transfers  electrons  to  the  transmission  line  as  it  propagates. 
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Fig.  10  Effective  velocity  matching  on  the  optical  input  response 
At  any  x ,  the  photon  density  ( per  unit  area )  remaining  in  the  light  wave  is 

F(x)  =  F,-e-r*  (23) 

where 

E0 
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and  is  moving  with  velocity  v0 .  Then  the  total  charge  per  unit  area  on  the  transmission 
line  is 

where 

Qo  =  mFo 

and  is  moving  with  velocity  ve .  The  current  produced  by  the  moving  charge 
is 

=  (25) 
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where  ac  represents  the  attenuation  of  the  electrical  wave  on  the  transmission  line.  From 
the  relationship 
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By  using  the  Fourier  transform,  we  can  convert  time  domain  into  frequency  domain,  then 
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We  can  assume  the  input  current  is 

/(»,/)= 

Then,  the  modulation  transfer  function  is  expressed 

HiA-lM. 


Ya-e~a’le 


(Ta-  ae)-  jlxf 


V  —  V  V 

v  o  y  e  y  o 


V  V  V 
o  y e  J  ye 


V  V 

-jltfJtSJL. 

v  -v„ 


28 
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It  is  noted  that  x  is  introduced  as  a  measure  of  velocity  mismatch  and  if  there  is  no 
mismatch  (v0  =  ve) ,  then  x  is  undefined  .  In  that  case  (16)  is  no  longer  applicable  and  the 

current  can  no  longer  be  defined ,  i.e.  if  an  ideal  delta  function  of  shape  is  introduced  and 
there  is  no  pulse  spreading,  the  current  response  is  infinite  which  is  the  limit  approached 
by  (17).  Equation  (18)  defines  the  impulse  response  of  the  system  as 


m=W)=e  ' '  '  7^7 

It  is  noted  that  ae  and  ve  are  frequency  dependent  and  therefore  the  pulse  response  will 
degrade  as  inputs  with  higher  and  higher  frequency  content  are  considered.  To  show  these 
effects  more  clearly,  it  is  more  appropriate  to  consider  the  transfer  function 
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Therefore  the  bandwidth  of  the  detector  is  determined  by  the  amplitude  and  phase  of 
H(f)  given  by 
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The  frequency  dependence  of  //(/ )  is  established  through  that  of  ae  which  is  decided  by 

(4),  (5),  and  (6),  using  (7)  and  (8).  The  frequency  dependence  of  the  group  velocity  ( or 
energy  velocity  )  ve  is  given  by 
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with  seff  ( co )  as  determined  by  (6) 


The  detector  transfer  function  is  plotted  in  Fig.l  1.  the  3  dB  frequency  of  the  detector  is 
given  by 


and  is  clearly  reduced  by  the  frequency  dependence  of  both  ae  and  a0 .  It  is  interesting  to 
note  that  the  bandwidth  is  independent  of  the  length  of  the  detector  which  results  because 
of  the  matched  load  configuration. 


TRAVELING  WAVE  DETECTOR  DYNAMIC  RESPONSE 


Fig.l  1  Traveling  Wave  Detector  Dynamic  Response 
8.2.2  Parasitic  Time  Constant 

In  the  TW  HFET  detector,  the  resistance  from  the  absorbing  channel  to  the 
contact  ( which  includes  the  doped  channel  resistance,  the  n+  channel  resistance  and  the 
metal  contact  resistance )  combine  with  the  capacitance  of  the  doped  channel  to  gate 
contact  to  produce  a  parasitic  RC  time  constant.  This  time  constant  constitutes  delay  in 
adding  charge  to  the  transmission  line  which  results  in  an  effective  velocity  mismatch 
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imposed  before  the  onset  of  transmission  line  propagation.  The  3dB  response  which  is 
also  shown  in  Fig.  10  would  be  modified  to 


where  o)m  is  given  by  (33)  and 
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9.  FABRICATION  AND  PHOTOMASK  DESIGN  AND  GENERATION 


The  fabrication  sequence  was  redesigned  to  implement  traveling  wave  devices. 
Several  changes  were  made  as  follows 

•  The  Si  implant  will  be  used  for  the  laser  and  detector  current  steering 
function. 

It  will  also  be  used  to  create  the  proper  threshold  voltage  for  depletion  transistors.  This  is 
a  light  field  mask. 

•  Waveguides  must  be  formed  and  the  waveguide  core  must  be  disordered  to 
shift  its  bandgap  to  higher  energy  to  produce  low  loss.  To  achieve  this,  the  Si  implant  that 
forms  the  source  and  drain  regions  is  now  performed  as  a  blanket  implant  over  the  entire 
wafer.  Where  the  Si  is  required  in  the  waveguide  for  disorder,  it  remains  and  where  it  is 
required  as  the  source  and  drain  of  active  devices  it  remains.  Everywhere  else  it  is 
removed  after  the  RTA  step.  The  RTA  step  is  performed  both  to  activate  source/drain 
regions  and  to  disorder  the  waveguide.  The  implant  is  performed  everywhere  except 
where  blocked  by  W  or  resist. 

•  AlxOy  is  used  for  isolation ,  passivation  and  to  form  the  lower  bragg  reflector 
for  the  optical  devices.  However  it  is  not  used  to  form  the  current  funneling  layer  as  in 
most  approaches  today.  Placing  AlAs  that  close  to  the  interface  would  compromise  the 
operation  of  the  FET.  We  see  no  possibility  to  use  the  oxide  to  form  a  high  performance 
FET.  The  HFET  or  PHEMT  is  already  the  ideal  FET.  Further,  the  concept  of  oxidizing 
large  lateral  distances  to  produce  a  small  aperture  is  not  manufacturable.  That  is  why  the 
Si  implant  step  is  chosen  above.  Since  we  are  using  oxidation  of  the  bottom  mirror,  then 
an  etch  step  referred  to  as  trench  etch  is  added  which  etches  through  the  back  mirror  prior 
to  oxidation.  Thus  all  devices  are  isolated  by  the  trench  etch  and  after  oxidation  all  the 
etched  vertical  surfaces  are  passivated  by  oxide. 

It  is  noted  that  the  oxide  is  used  primarily  to  produce  the  lower  mirror.  Since  the 
oxide  rate  is  quite  sensitive  to  the  A1  concentration  then  the  AI0.7  layers  will  oxidize  much 
more  slowly  than  the  A1  layers.  Therefore  when  20-30  pm  of  AlAs  is  oxidized  in  the 
bragg  mirror,  there  may  be  only  2  pm  of  oxidation  of  AI0.7.  Since  this  oxidation  occurs 
under  the  source  and  drain  regions  and  along  the  waveguide  walls  then  it  serves 
additional  purposes  which  are  1)  it  reduces  the  capacitances  of  the  sources  and  drain 
regions  to  produce  higher  speed  FETs  2)  it  provides  a  low  index  lateral  cladding  on  the 
waveguide  to  allow  very  tight  waveguide  bends.  Just  how  far  we  can  push  this  concept  of 
lateral  oxidation  remains  to  be  determined  experimentally. 

We  have  devised  a  unique  approach  to  use  the  selective  oxidation  of  A1 
containing  compounds  over  GaAs  to  produce  our  structures  which  follows  the  same 
principles  as  the  Si  LOCOS  technology.  The  oxidation  is  done  with  GaAs  as  a  mask,  with 
W  as  a  mask  and  with  Si02  as  a  mask. 

With  these  modifications  the  process  flow  and  masking  steps  are  as  follows: 

1.  Pattern  mask  1  to  etch  alignment  marks.  Deposit  protective  oxide  over  apertures 

2. Pattem  mask  2  for  first  Si  implant  and  perform  implant 
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3.  Pattern  for  lift  off  of  W.  Etch  Si02.  Sputter  W  and  lift-off 

4. Pattem  W  gates,  etch  in  SF6 

5.  Blanket  ion  implant.  RTA  for  disorder  and  activation. 

6. Pattem  Source/Drain  implants  to  protect  N+ junctions  and  waveguides.  Etch  to  top  of 
AI0.7  layers. 

7. Pattem  for  trench  etch.  Etch  trench  and  then  oxidize 

8. Pattem  to  protect  waveguide,  etch  to  etch  stop  in  S/D  areas 

9. Pattem  for  subcollector  etch,  etch  oxide  and  semiconductor  with  RIE 

10. Pattem  for  N  lift-off,  evaporate  and  lift-off 

1 1  .Pattern  for  collector  etch  ,  etch  oxide  and  semiconductor  with  RIE 

12. Pattem  for  P  lift-off,  evaporate  and  lift-off 

13.  Apply  polyimide  isolation 

14. Pattem  for  polyimide  holes 

15.  Lift-off  of  final  Au  bond  pads 

16. Deposit  first  stack  layer.  Pattern  grating,  and  RIE  etch  grating  pattern.  Finish  stack 
deposition 

17. Etch  stack  for  final  Au  bond  pads 

The  total  mask  count  is  15  or  14  if  we  choose  to  delete  the  pattern  in  step  8.  This 
becomes  a  question  of  how  the  waveguide  losses  are  affected  by  forcing  the  mode  close 
to  the  top  mirror.  It  probably  can  be  deleted. 

The  mask  set  was  designed  for  maximum  debugging  capability.  Several  chip  sites 
were  developed.  One  of  these  was  shortened  to  7  levels  to  produce  only  FET  devices. 
Another  chip  site  was  shortened  to  only  8  levels  to  produce  only  VCSELs.  These  are  top 
emitting  rectangular  devices  of  various  sizes  and  aspect  ratios.  Another  chip  site  was 
devoted  to  process  diagnostics  including  mirror  testing  areas  and  also  etch  depth  testers 
for  all  of  the  various  etches  involved  in  the  process.  The  purpose  of  the  shortened 
fabrication  chip  sites  was  to  provide  a  means  to  obtain  feedback  on  laser  or  FET  devices 
independently  with  fewer  man  hours  than  required  for  the  grating  laser.  These  modules 
need  to  work  first  to  verify  the  integrity  of  the  approach .  When  these  have  proven 
themselves  we  will  implement  the  full  process  to  make  the  grating  laser. 

The  laser/detector  requires  the  grating  level  which  is  a  non-standard  item.  We 
have  chosen  the  phase  mask  approach  to  achieve  this.  The  phase  mask  can  be  made  by 
holography  or  by  e-beam  printing  of  each  desired  feature  and  then  transferring  these 
patterns  to  a  quartz  plate.  The  plate  is  used  as  a  photomask  and  exposed  with  248  nm 
light.  We  worked  with  two  companies  to  make  the  phase  mask,  one  using  e-beam  charged 
$10.5K  and  the  other  using  holography  (Denmark)  charged  S5.5K.  It  seems  the 
holography  approach  is  correct  for  initial  testing  because  the  grating  pitches  on  all  devices 
are  the  same  and  the  devices  have  been  arranged  to  be  in  parallel. 

A  major  activity  was  invested  in  the  development  of  the  photomask  set  involving 
more  than  12  man  months  of  effort.  We  were  bridging  the  gap  between  the  earlier  vertical 
cavity  process  flow  that  only  begun  to  yield  results  at  AT&T.  Completely  reworking  this 
procedure  to  incorporate  waveguides ,  traveling  wave  structures,  a  new  oxidation 
technique,  disorder  etc.  turned  out  to  be  an  enormous  task.  The  masks  were  produced 
with  the  best  simulation  support  we  could  provide.  These  masks  contain  a  very  wide 
variety  of  circuits  and  devices.  The  list  of  circuits  and  devices  includes  among  others 
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ring  oscillators  to  determine  transistor  speed 

ring  oscillators  to  determine  switch  laser  speed 

12  traveling  wave  detectors  (geometry  and  output  variations) 

12  traveling  wave  lasers  (geometry  and  output  variations) 

a  traveling  wave  DOES  detector  receiver 

temperature  tuned  lasers  and  detectors  using  FET  control 

waveguide  combined  lasers  with  temp,  tuning  to  produce  1 00GHz  optical  RF 

a  directional  coupler  using  the  HFET  modulator  with  dual  waveguide  I/O 

optical  waveguide  amplifiers 

HFET  CCD  arrays  with  optical  outputs 

numerous  test  devices  for  the  DOES,  HFET  laser  and  HFET  transistors 
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10.  GROWTH  AND  TESTING 


The  program  required  that  no  more  than  20%  of  the  effort  be  devoted  to 
simulation  but  in  our  case  this  turned  out  to  be  impossible.  In  April  1994  the  PI  launched 
an  effort  at  UConn  to  install  a  state  of  the  art  clean  room  and  growth  facility  which 
included  acquisition  of  space  and  funds  and  the  design  and  implementation.  Due  to  the 
bidding  process  and  uncountable  unforeseen  delays,  construction  did  not  finally  get 
under  way  until  June  1996  and  the  facility  was  occupied  in  Nov.  1996.  With  the 
equipment  installation  and  required  qualification  procedures  we  ran  out  of  time  to 
produce  working  devices.  Another  major  problem  was  the  supply  of  MBE  material.  A 
Varian  GEN  II  system  has  been  installed  in  the  new  growth  facility  at  UConn  but  it  was 
not  on  line  at  that  time  to  produce  the  wafers  for  this  program.  Thus  an  arrangement  was 
made  with  Penn  State  (PSU)  to  supply  wafers.  The  PI  had  been  awarded  an  SBIR  (Rome 
Labs)  contract  with  a  company  in  MD  to  implement  neurodes  for  an  optoelectronic 
processor  and  PSU  was  to  receive  $28K  per  year  for  wafers.  The  PI  was  also  trying  to 
raise  venture  capital  to  continue  this  work.  At  the  time  of  the  Phase  II  award  (Sept  1996), 
a  conflict  arose  as  to  the  ownership  of  IP  between  the  small  company  and  the  venture 
people.  Due  to  an  unfortunate  turn  of  events  the  small  company  took  the  money 
elsewhere  and  this  source  of  funding  disappeared.  Thus  PSU  was  unwilling  to  supply 
wafers  without  funds.  By  the  end  of  the  contract  (March  1997)  we  were  still  trying  to 
obtain  wafers.  The  MBE  reactor  at  UConn  is  still  being  qualified  but  will  require  an 
infusion  of  funds  before  serious  MBE  growth  can  begin.  In  order  to  perform  the  testing 
the  PI  acquired  an  HP  design  and  test  station  at  a  total  cost  of  $25  OK  that  will  perform  on 
wafer  measurements  up  to  50GHz.  About  $25K  of  contract  money  was  used  to  start  the 
lease  for  this  equipment  and  the  PI  is  now  trying  to  raise  the  funds  to  continue  this  lease 
within  the  university  and  by  finding  new  funding  sources.  Therefore  during  this  one  year 
program,  the  PI  did  all  of  the  fundamental  design  work  to  support  this  new  concept,  and 
he  assembled  all  of  the  necessary  tools  and  laboratories  to  continue  this  work. 

At  the  present  time ,  funding  is  a  major  problem.  A  new  SBIR  has  been  started 
(Phase  I)  to  demonstrate  true  time  delay  circuits  for  phased  array  systems.  This  contract  is 
with  the  NavAir  department  of  the  DoD.  Any  results  arising  from  this  SBIR  will  be  made 
available  to  the  airforce  electro-optics  directorate  on  a  continuous  and  timely  basis.  The 
Phase  I  work  will  result  in  a  small  fabrication  activity  to  build  a  directional  coupler. 
Hopefully  the  Phase  II  will  enable  a  substantial  effort.  Also  recently  an  STTR  was 
awarded  to  demonstrate  a  CCD  with  infra  red  detection  using  the  QWIP  from  the 
BMDO.  This  may  provide  another  source  if  it  can  go  to  a  Phase  II. 

The  attempt  to  raise  venture  capital  was  not  really  successful.  They  agreed  to 
supply  $360K  but  they  wanted  another  investor.  We  could  not  find  another  investor.  An 
ATP  was  applied  for  but  the  outcome  is  unknown.  NSF  proposals  have  been  submitted. 
More  SBIR  proposals  will  be  submitted.  Visits  have  been  made  to  DARPA  to  seek 
support  but  requests  for  funding  have  not  been  successful.  There  is  some  potential  interest 
from  Celeritek,  a  small  company  in  Ca.  which  manufactures  GaAs  circuits  for  wireless 
and  DBS  systems  (this  company  was  founded  with  help  from  Sutter  Hill  Ventures  and 
were  briefed  on  the  ICT).  They  have  determined  that  an  optical  link  from  the  satellite 
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receiver  to  the  set  top  would  be  an  exciting  product  and  may  be  interested  in  funding 
(such  a  link  would  allow  an  IF  frequency  of  4GHz  instead  of  70MHz  currently  and  could 
dramatically  increase  data  flow).  Only  a  monolithic  chip  could  achieve  a  low  enough  cost 
to  make  this  concept  viable.  There  is  a  possibility  that  Celeritek  may  provide  the  funds 
for  a  prototype  if  they  can  accept  the  risk.  Their  decision  should  be  known  within  two 
months. 
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ACRONYMS 


E/O,  O/E 

Optical  to  Electrical,  Electrical  to  Optical 

BER 

Bit  Error  Rate 

ATM 

Asynchronous  Transfer  Mode 

ICT 

Inversion  Channel  Technology 

MBE 

Molecule  Beam  Epitaxy 

PHEMT 

Pseudomorphic  High  Electron  Mobility  Transistor 

III-V 

Semiconductor  compounds  using  group  III  and  Group  V  elements 

SI  GaAs 

Semi-Insulting  GaAs 

IC’s 

Integrated  Circuits 

FET 

Field  Effect  Transistor 

DOES 

Double  Heterostructure  Optoelectronic  Switch 

HFET 

Heterostructure  Field  Effect  Transistor 

CAD 

Computer  Aid  Design 

UConn 

University  of  Connecticut 

SCH 

Separate  Confinement  Heterostructure 

Q  factor 

Quality  factor 

GC  VCSEL 

Grating  Coupled  Vertical  Cavity  Surface  Emitting  Laser 

VCSEL 

Vertical  Cavity  Surface  Emitting  Laser 

DFB 

Distributed  Feedback 

DBR 

Distributed  Bragg  Reflection 

PC 

Personal  Computer 

VC 

Vertical  Cavity 
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TMM 

DC/AC 

RIN 

RTA 

RC 

RCE 

TWPD 

WDM 

TW 

R 

G 

LOCOS 

S/D 

RIE 

PSU 

SBIR 

CCD 

QWIP 

DBS 

IF 


Transmission  Matrix  Method 
Direct  Current/ Alternating  Current 
Relative  Intensity  Noise 
Rapid  Thermal  Anneal 
Resistance  Capacity 
Resonant  Cavity  Enhanced 
Traveling  Wave  Photon  Detector 
Wavelength  Division  Multiplexing 
Traveling  Wave 
series  Resistance  per  unit  length 
shunt  conductance  per  unit  length 
Local  Oxidation  of  Silicon 
Source/Drain 
Reactive  Ion  Etching 
Penn  State  University 
Small  Business  Innovation  Research 
Charge  Coupled  Device 
Quantum  Well  Infrared  Detector 
Direct  Broadcast  System 
Intermediate  Frequency 
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LIST  OF  SYMBOLS 


Gb/s  Giga  bit  per  second 

ni ,  n3  the  index  of  the  cladding  region 

ri2  the  index  of  the  core  region 

LPb,  Lpt  the  penetration  depths  of  the  bottom  and  top  mirrors 

An  the  index  difference  between  the  two  layers  in  a  pair  of  the  mirror 

LMt,b  is  the  total  front  and  back  mirror  thicknesses 

k  the  coupling  constant 

Lx  the  quantum  well  width 

jje  the  electrical  confinement  factor 

J ,  JTH  the  current  and  theoreld  current  density 

rp  the  effective  photon  lifetime  for  the  vertical  cavity  device 

avc  the  vertical  cavity  loss  due  to  the  diffraction  grating 

Pouc  power  output  from  the  laser 

tf  the  power  transmissivity  of  the  output  port 

avc  the  primary  optical  loss  in  the  vertical  cavity 

ijd,ff  the  diffraction  efficiency  due  to  grating 

r|  i  ?  r|2  the  impedances  of  the  guide  and  the  incident  medium 

Ay  the  field  intensity  at  the  edge  of  the  guide  of  the  z  propagating  wave 

C0  the  field  intensity  of  the  incident  wave 

(Xg  the  growth  of  the  guided  wave  due  to  diffraction 

deff  the  effective  width  of  the  guide  which  includes  evanescent  cladding  penetration. 

ax  the  loss  of  wave  propagating  along  x  direction  (vertical) 

r  the  confinement  factor 

Fg  the  waveguide  photon  density 

t'p  the  effective  photon  lifetime  in  the  cavity 

W  the  width  of  waveguide, 

L  the  length  of  waveguide 
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rie  the  electrical  efficiency  of  laser 

Z0.  the  characteristic  impedance 

X  the  wavelength 

0)  the  frequency 

ws  ,  wg  the  source  and  gate  widths, 

q  the  permittivity 

cr.  ,  g„  the  metal  conductivities 

K.(k , ),  AT(Ar,  ),&(&)  the  complete  elliptic  integrals  of  the  first  kind 

the  effective  index 

E0  the  delta  function  input  impulse  of  energy 
a  the  absorption  coefficient  of  the  active  layer 

r  the  optical  confinement  factor 

v0  the  optical  velocity 

ve  the  electrical  energy  (or  group)  velocity, 

/  the  detector  length 

t  the  time  as  measured  from  the  time  of  incidence  of  energy  to  the  detector 
77,  the  internal  quantum  efficiency 

q  the  electronic  charge 

hv  the  photon  energy 

ac  the  attenuation  of  the  electrical  wave  on  the  transmission  line. 

H[f)  the  transfer  function 
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12.  Appendices 


Appendix  A.  Determination  of  Diffraction  Efficiency  for  a  second  order  Corrugated 
Waveguide 

Determination  of  Diffraction  Efficiency  for  a  Second 
Order  Corrugated  Waveguide 

G.  W.  Taylor  and  C.  Kwan 
Electrical  and  Systems  Engineering 
University  of  Connecticut 
Storrs,  CT  06269 

Abstract 

The  diffraction  from  a  second  order  corrugated  waveguide  is  analyzed  to 
determine  its  efficiency  as  an  element  to  redirect  a  guided  wave.  Using  a  rect¬ 
angular  grating,  the  efficiency  is  evaluated  in  an  approximate  closed  form.  The 
dependence  upon  grating  thickness,  waveguide  thickness  and  wavelength  are  in¬ 
vestigated  and  the  results  indicate  usefulness  as  a  simulation  tool  for  photonic 
integrated  circuits. 
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1  Introduction 


Gratings  combined  with  waveguides  have  been  important  structures  for  realizing 
many  useful  components  for  integrated  optics  [1]  [2].  These  include  co-directional 
and  contra-directional  waveguide  couplers  for  optical  switches  [3]  and  DFB  (Dis¬ 
tributed  feedback)  or  DBR  (Distributed  Bragg  Reflector)  laser  structures  [4].  In  the 
application  to  lasers,  first  and  third  order  gratings  are  normally  applied  to  make  clas¬ 
sic  DFB  or  DBR  devices  since  these  gratings  produce  diffracted  components  only  in 
the  plane  of  the  guide.  However  second  order  gratings  have  been  considered  as  a 
means  to  discriminate  between  the  two  symmetrical  modes  produced  by  a  DFB  laser, 
i.e.  the  second  order  grating  produces  components  in  forward,  backward  and  surface 
normal  directions  and  the  lasing  action  will  naturally  select  the  mode  which  favors 
the  least  surface  normal  radiation  loss  [5]  [6].  Thus  the  wavelength  is  selected  and 
stabilized.  In  the  past  few  years,  second  order  gratings  have  been  used  as  the  output 
coupling  device  for  DBR  lasers  to  produce  surface-emitting  devices  either  individually 
(as  a  broad  area  laser)  or  as  arrays  to  generate  significant  power  [7]  [8].  Another  use 
of  the  second  order  grating  [9]  is  to  couple  from  a  passive  waveguide  into  an  active 
optoelectronic  component  such  as  a  detector  or  a  modulator  or  of  coupling  the  light 
into  a  passive  waveguide  from  a  laser.  The  most  significant  parameter  characterizing 
such  grating  structures  is  the  loss  parameter  a  which  describes  the  fraction  of  power 
propagating  in  the  guide  which  can  be  coupled  to  the  surface  normal  direction.  In 
general,  this  parameter  has  been  determined  from  numerical  simulations.  It  is  the 
purpose  of  this  paper  to  investigate  the  calculation  of  a  and  its  dependence  on  the 
waveguide  and  grating  parameters  using  quasi-analytical  techniques. 


2  Second  order  grating  calculation 


The  original  calculations  which  formalized  Maxwell’s  equations  to  address  this  prob¬ 
lem  were  published  by  Streifer  et  al.  [10].  The  approach  was  generalized  to  any  order 
and  to  any  shape  of  grating.  A  significant  result  of  their  work  was  a  numerical  al¬ 
gorithm  which  could  be  used  to  determine  the  behavior  of  the  grating.  Subsequent 
calculations  by  Streifer  [11]  and  others  [5]  have  focussed  on  the  analysis  of  the  second 
order  grating  in  the  context  of  the  DFB  laser  in  which  coupled  mode  analysis  was 
employed  to  determine  the  behavior  of  the  grating  in  the  presence  of  symmetrical  for¬ 
ward  and  backward  travelling  coupled  modes.  Here  also  the  analysis  yields  numerical 
solutions.  Our  interest  here  is  to  re-examine  the  basic  diffraction  for  a  wave  inci¬ 
dent  upon  a  second  order  grating  and  to  determine  an  analytic  form  for  the  grating 
coupling  coefficient  in  terms  of  the  structural  parameters.  By  representing,  the  first 
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order  backward  travelling  wave  as  a  reduced  version  of  the  forward  wave,  we  are  able 
to  consider  both  weak  and  strong  coupling  situations.  This  enables  a  simplification 
and  thus  an  analytic  result  is  possible.  The  result  we  obtain  will  be  useful  in  evalu¬ 
ating  and  designing  waveguide  couplers  as  inputs  or  outputs  for  lasers  and  detectors 
into  waveguides  without  resorting  to  extensive  numerical  methods.  These  simulated 
results  are  ideally  suited  as  the  basis  of  CAD  tools  for  integrated  optic  design. 


The  problem  we  address  is  illustrated  in  Fig.  1.  A  slab  waveguide  of 
thickness  d  is  corrugated  with  a  grating  of  thickness  t.  The  grating  begins  at  z  =  0 
and  extends  continuously  in  the  positive  z  direction.  The  x  direction  is  taken  as 
positive  downwards.  The  slab  waveguide  is  characterized  by  the  indices  nx  above,  n2 
within,  and  n3  below  the  waveguide  respectively.  The  grating  section  has  a  refractive 
index  with  a  periodic  variation  and  is  described  by  the  Fourier  series 

OO 

n2g{x,z)  =  n2(x)+  Bq{  x)el2nqz/A  (1) 

qc=—oo  &  q^O 

where  n  is  the  average  index  over  one  grating  period,  Bq  are  the  Fourier  coefficients, 
and  A  is  the  grating  period  as  illustrated  in  Fig.  1. 


The  TE  modes  in  the  waveguide  are  described  by  the  two  dimensional 
wave  equation 

8PEy  SPEy  l2  O,  ,ns 

~dx?  +  ~df+  °n  (*’ z)  y  =  (2) 

where  A;0  =  27t/Ao  is  the  wave  number,  and  A0  is  the  free  space  wavelength.  The  index 
is  expressed  by  region  as 


n2(x,z)  =  ^ 


nl  x  <  0 

n2(x,z)  0  <x<t 

n|  t  <  x  <d 

rvi  d  <  x 


(3) 


In  the  standard  approach  to  this  problem  the  mode  field  Ey  is  represented  as  shown 
by  Streifer  [10],  by  an  infinite  sum  of  partial  waves  in  the  form 

OO 

Ey(x,z)=  Em(x,z)exp(ipmz)  (4) 

m=—oo 


where  the  propagation  constants  (3m  are  expressed 


Pm  =  Po  ~ 


2tt  m 

IT 


(5) 


and  we  note  that  (5)  implies  (as  shown  later)  that  m  takes  positive  values  for  diffracted 
waves.  If  no  grating  were  present  then  all  Em(x,z)  with  would  be  zero,  and 
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we  would  have  E0(x)  and  /3q  as  the  mode  pattern  and  propagation  constant  of  the 
unperturbed  waveguide.  Therefore  in  our  analysis  we  assume  that  a  wave  Ein 

Ein  =  Mx)e~i0o=  (6) 


is  propagating  in  the  waveguide  for  z  <  0  where  A0(x)  is  an  arbitrary  amplitude.  We 
assume  a  TE  mode  in  the  guide  and  therefore  A0(x)  is  determined  by 


Aq(x) 


'  AyeSx  x  <0 

<  Ay  cos  kx  +  By  sin  kx  0  <  x  <  d 

(Ay  cos  kcL  +  By  sin  Kd)e~Ax~d)  d  <  x 


(7) 


where 


k  =  n^ko  sin  9b, 

6  =  (n|  —  nf)  k2Q  -  nfal  sin2  9b 

7  =  («!  ~  nl)  ko  ~  ™2ko  sin2  °b 
Ay  is  the  value  of  A)(0)> 


1/2 

1/2 


and  9b  is  the  bounce  angle. 


From  continuity  of  Hs  at  x  =  0,  we  also  obtain 


Ay  _  K 


Bv 


(8) 


For  the  TE  mode,  the  graphical  solution  of  the  eigenvalue  equation,  ta nnd  = 
yields  the  eigenvalue  k.  From  the  eigenvalues  we  also  obtain  the  bounce  angles  for 
TE  modes  as 

9b  =  sm  -  (9) 

K 

and  from  basic  waveguide  theory  we  know  that  the  propagation  constant  fio  is  related 
to  the  bounce  angle  9b  by 

p0  =  n2k0  cos  9b  (10) 


Our  analysis  will  be  confined  to  the  problem  of  the  second  order  grating 
imposed  on  the  waveguide.  This  grating  is  of  particular  importance  to  applications 
which  require  the  redirection  of  optical  waves  through  90°  such  as  in  coupling  from 
a  waveguide  to  a  detector.  Second  order  grating  diffraction  principles  are  outlined  in 
many  texts  (e.g.[12j)  and,  as  illustrated  in  Fig.  2a,  the  light  is  diffracted  by  an  angle 
^ )  according  to  the  Bragg  relation 


—cosip  =  sin9  =  m  —  1 


(11) 
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where  ip  =  9  +  7r/2  is  the  angle  of  propagation  of  the  diffracted  wave  with  respect 
to  the  direction  of  propagation  of  the  incoming  wave.  For  real  8,  only  the  values  of 
m  —  0,  1,  and  2  are  allowed  corresponding  to  components  in  the  forward  direction 
(m  =  0),  the  reverse  direction  (m  =  2)  and  the  direction  approximately  normal  to 
the  waveguide  (m  —  1).  For  9  =  0,  the  wave  propagates  normal  to  the  guide  and 
for  small  positive  or  negative  values  it  propagates  to  the  right  or  left  of  the  normal 
respectively  of  the  guide  in  Fig.  2a.  The  component  approximately  normal  to  the 
waveguide  is  comprised  of  upward  and  downward  travelling  waves  since  both  will  be 
allowed  by  the  Bragg  relation.  Their  relative  magnitudes  will  be  determined  by  the 
indices  involved  in  a  particular  case  as  illustrated  by  example  later.  Therefore  in 
consideration  of  (4)  only  three  terms  in  the  summation  need  to  be  considered.  The 
diffraction  orders  can  be  shown  with  the  vector  diagram  as  illustrated  in  Fig.  2b. 

We  first  consider  the  rectangular  grating  because  of  its  mathematical  sim¬ 
plicity.  In  the  appendix,  the  sinusoidal  grating  is  considered 


3  Rectangular  grating 


The  rectangular  grating  has  been  discussed 
the  coefficients  in  (1)  are  expressed 

("l  -  ) 


B,= 


0 


7r  q 


A 


in  detail  by  Streifer  [10]. 


0  <  x  <  t 

x  <  0,  t  <  x,  gf^O 


For  this  case 


(12) 


where  w  is  the  grating  tooth  width  (see  Fig.  4).  The  index  in  the  grating  region 
corresponding  to  (3)  is  expressed 


n2g(x,z,  q ) 


n2 


M  -  tg) 

qn 


sin( 


7T  qw 

~A~ 


where  the  average  n  is  given  by 

ri  =  {wnl  +  (A  —  u>)n|}/ A  ,  0  <  x  <  t 


(13) 

(14) 


The  remaining  indices  axe  as  in  (3).  The  first  step  in  the  solution  is  to  determine 
Ey(x,  z)  in  the  region  of  the  grating  itself.  However,  we  first  note  that  the  coefficients 
q  —  — 1  and  +1  determine  the  coefficients  for  the  waves  with  m  =  0  and  2  respectively. 
Separate  wave  equations  are  obtained  for  each  value  of  m  by  substituting  (4)  into  (2) 
and  collecting  the  coefficients  for  each  value  of  et/3m~.  We  obtain 

-  $E0  +  k20n2E0  +  2z/?0^  +  k20  (Bt  El  +  B2  E2)  =  0,  m  =  0  (15) 

oxz  oz 
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(16) 


o2  T? 

^  -  fiEi  +  kin2 Ei  +  2 +  k2  (B-i  E0  +  Br  E2)  =  0,  m  =  1 

^  -  (32E2  +  k2n2E2  +  +  k2  (B-i  Ei  +  B_2  E0)  =  0,  m  =  2  (17) 

where  we  represent  S_i  =  Bx  =  B  =  -(^~n^w(j)  as  determined  by  (13)  and, 
as  is  normally  done,  the  terms  in  are  neglected.  It  is  noted  that  these  equations 
more  closely  resemble  those  of  Kazarinov  and  Henry  [5]  rather  than  Streifer  [10]  since 
Em(x,z )  in  (4)  has  been  given  an  explicit  z  dependence  which  results  in  the  extra 
term  of  dEi/dz  in  (16).  This  term  will  be  important  in  considering  diffraction  at 
angles  of  'ip  other  than  7t/2  (i.  e.  for  non-zero  values  of  6). 

These  are  three  coupled  equations  in  the  unknowns  Eq ,  Ei  and  E2  and 
it  has  been  shown  previously  by  Streifer  and  others  that  the  only  way  of  obtaining 
solutions  is  through  a  numerical  perturbation.  For  the  perturbation,  the  term  in  B 
in  each  equation  is  initially  neglected  to  obtain  a  first  attempt  solution.  These  results 
are  then  used  to  determine  the  B  terms  and  the  procedure  repeated  until  convergence 
is  obtained.  Boundary  conditions  are  also  applied  .  For  each  of  the  solutions  Eq,  Ei 
and  E2,  the  continuity  of  E  and  dE/dx  must  be  applied  at  the  interfaces  x  =  0,  x  =  t 
and  x  —  d.  It  is  thus  required  to  obtain  the  solution  of  18  equations  in  18  unknowns 
with  one  extra  relation  for  input  power  conservation. 


3.1  Ideal  Second  Order  Grating 

To  circumvent  this  lengthy  numerical  procedure,  a  closed  form  result  may  be  obtained 
with  some  reasonable  approximations.  In  Fig.  3  we  sketch  the  nature  of  the  solution. 
The  wave  E{n  enters  the  grating  from  the  left  and  decays  with  the  constant  a  as  it 
travels  through  the  grating,  giving  rise  to  the  two  guided  waves  E0e~a::  and  E2eaz 
and  a  radiating  wave  E\.  It  is  clear  that  the  wave  Ex  arises  from  the  losses  that  are 
incurred  by  Eq  as  it  propagates  through  the  grating.  Although  Ei  represents  a  wave 
propagating  in  the  positive  and  negative  x  directions  it  will  also  ,  in  general,  have  a  z 
dependence  resulting  in  a  finite  dEi/dz  in  (16).  However  for  the  second  order  grating 
it  is  true  in  general  that  (3q  —  —j32  and  then  it  can  be  shown  from  (5)  that 

Pi  =  0  with  po  =  27t/A  (18) 

Therefore  the  terms  P\E\  and  2i(3idEi/dx  disappear  in  (16).  Equations  (10)  and  (18) 
state  that  for  a  perfect  second  order  grating  the  grating  pitch  is  fairly  close  to  the 
wavelength  and  is  dependent  on  the  bounce  angle  6b  and  therefore  on  the  thickness  of 
the  guide  d,  through  (9).  For  a  two  dimensional  or  rib  waveguide  /3q  will  also  depend 
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_ S2B 


on  the  lateral  geometries.  Therefore  it  is  necessary  to  know  the  exact  waveguide 
shape  before  designing  the  grating  for  second  order  diffraction. 

The  grating  produces  a  first  order  reflected  wave  E2ej/3°z  travelling  in  the 
— z  direction,  as  Figure  3  shows,  which  is  a  diminished  version  of  Eq  in  amplitude. 
Since  it  is  only  the  amplitude  is  required  in  (16),  then  E2  can  be  represented  from 
transmission  matrix  theory  as 


E2  =  Eonmeff,  (rg  =  rxmeff ) 

where  rx  is  the  reflectivity  of  one  grating  period,  rg  is  the  total  reflectivity  of  the 
grating,  and  meff  —  is  the  effective  number  of  mirror  periods  seen  by  the 

incident  field.  For  a  strictly  first  order  grating  it  can  be  shown  that  rg  is  a  real 
number  and  approaches  unity  for  an  infinite  grating,  i.e. 

r  i 

5  i  +  K/^)2m 

where  m  is  the  number  of  grating  periods.  Due  to  the  existence  of  the  second  order 
component,  rg  will  be  less  than  unity  but  we  can  still  express  the  z  propagating 
grating  field  as 

E0(x,  z)  =  E0e ~az  +  rgE0e~az  =  E0(  1  +  rg)  e~«z  (19a) 


An  important  approximation  we  make  regards  Eq  and  E2.  Since  Eq  is  a 
continuation  of  Ein  and  E2  is  a  reflection  of  Eq,  then  instead  of  finding  Eq  and  E2 
in  the  grating  and  matching  it  to  the  remainder  of  the  guide  (i.e.  finding  Eq  and  E2 
solutions  in  all  sections  of  the  guide),  we  represent  E0+  E2  as  one  solution  for  the 
entire  guide  in  the  form  of  Eq  as  given  by  (19a).  These  two  waves  are  assumed  to 
continue  to  propagate  in  the  grating  with  the  constants  /50  and  — /30,  where  (5q  is  the 
propagation  constant  of  the  input  wave.  Obviously,  the  actual  propagation  constant 
will  be  somewhere  between  that  for  a  waveguide  of  thickness  d  and  thickness  d  —  t 
and  as  the  depth  of  the  grating  t  increases  and  begins  to  approach  the  waveguide 
thickness  d,  this  assumption  no  longer  becomes  tenable.  This  issue  is  discussed  in  the 
context  of  the  results.  The  x  dependence  of  the  Eq  wave  given  by  (19a)  is  contained 
in  the  waveguide  mode  profile  given  by  (7).  It  will  be  noted  therefore  that  Aq(x) 
and  thus  Ay  as  it  is  used  below  represents  the  net  amplitude  of  a  wave  travelling 
in  the  corrugated  section  of  the  guide  which  allows  us  to  calculate  the  diffraction 
efficiency  accurately  without  specifying  the  relative  magnitude  of  Eq  and  E2.Thus 
(19a)  becomes 

Eq{x,z)  =  A0(x)  e~az  (19b) 
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With  these  assumptions  we  need  to  solve  only  (16)  rewritten  as 


$2E  _  _ 

+  kffiEi  +  h%  B  E0(x,  z)  =  0  (20) 

together  with  (19).  Prom  the  solution  of  (20)  it  is  then  possible  to  determine  a,  which 
is  a  real  constant  and  can  be  found  from  the  application  of  boundary  conditions.  This 
is  a  major  distinction  from  previous  work  in  which  a  was  evaluated  from  numerical 
solutions  of  Ei  and  E0.  Since  Ei(x,  z )  is  the  solution  of  interest  we  refer  to  it  hence¬ 
forth  as  E  and  determine  its  solutions  E\x,  Eg,  E2x,  E3 x  in  the  various  regions  of  the 
guide.  This  approximation  may  be  summarized  by  stating  that  outside  the  grating, 
the  waveguide  is  a  three  layer  waveguide  and  that  within  the  grating  we  maintain  a 
three-layer  description.  Equation  (20)  then  yields  the  complementary  solution 

Egc{x )  =  C2eik°x  +  Cze~ik°x  ,  kg  =  k0n  (21) 

and  the  particular  solution  (see  Appendix  A), 

k2  ~  8  8 

Eov  =  o  °  ,  o  B  Ave~az(cos  kx  - I —  sin  kx)  =  b(z)  (cos  kx  -\ —  sin  /car)  (22) 
K?  -  k2  K  K 

which  defines  b(z)  =  ^^2  B  Aye~c **  =  b.  Then  the  total  solution  for  Eg  is 

C 

Eg(x,  z)  =  C2elkgX  +  Cze~lkaX  +  6( cos  kx  H —  sin  kx)  (23) 

AC 

The  solutions  in  the  non-grating  sections  of  the  waveguide  are  determined  from  (16) 
and  (3)  to  be  (note  that  f3\  =  0) 

Elx(x,  z)  =  Cielkix  ,  ki  =  (k^n2  -  02)l/2  =  k0rii  (24) 

E2x(x,  z)  =  C4eik^d~x)  +  C5e-ik2^  ,  k2  =  (k2r%  -  02)1/2  =  kQn2  (25) 

E3x(x,z)  =  C&e~lk3{x~d)  ,  k3  =  (k^nl  -  0\)l/2  =  k0n3  (26) 

The  arbitrary  constants  in  (23)-(26)  are  determined  by  the  continuity  of  E(x)  and 
dE/dx  (continuity  of  Hz)  at  the  boundaries  x  =  0,  x  =  t  and  x  =  d.  Then  the 
following  relations  are  obtained 

C1  =  C2  +  Cz  +  b,  x  =  0  (27) 

ik\C\  =  ikgC2  —  ikgCz  +  8  -b  ,  x  —  0  (28) 

C2eikgt  +  Cze~ikgt  +  b(c os  Kt  +  -  sin  Kt)  =  +  C5e~ik2^  ,  x  =  t  (29) 

K 
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ikg(C2eik<)t  —  Cz&  tkst)  — b(n sin Kt  — 6 cos Kt)  — —ik2(C4elk2^d  0  —  C$e  lk2(d  t'))  ,  x  =  t 

(30) 

C4  +  C5  =  C6  ,  X  =  d  (31) 

— ik2C4  +  ik2C$  —  — ikzCe  ,  x  =  d  (32) 


This  set  of  linear  equations  is  solved  in  a  straight-forward  way  and  some 
of  the  important  results  are 


C2  =  1(1  +  r-)Ct  ~  +  J-A  c*  =  ;(1  -  r)c>  -  5(1  -  (») 


kn 


ika 


and 


where 


Ci  =  26 


A  +  iB 

c+Jd 


=  2br(kg), 


C6  =  26{ 


A  +  iB 
C~+iD 


Yi  +  Yt} 


(34) 

(35) 


6  6 

A  =  cos (kgt)  +  —  sin (kgt)  -  [cos (nt)  +  -  sin(/ct)] 
kg  k 

1  +  ^){cos[2A:2(d  -t)  +  kgt]  +  £  sin[2 k2(d  -t)  +  kgt]} 

2  K2  —  k>Z  K2  Kg 

-  £){«*[»>(<*  -*)-  M  -  r  (<*  - 1)  -  M> 

2  rC  2  /C3  A& 

+(cosKt+  —  sin  /ct)  ( f 2  )  cos[2fc2  (d  —  t)] 

k  k2  —  k2 

+(j-  cos  Kt  —  sin  Kt) ( f 2  +  f 3 )  sin[2fc2(d  -  <)] 

k2  k2  k2  —  kz 

Jc  S  6  k 

B  =  -T1  sin(fcst)  +  —  cos  (kgt)  -  [—  cos  (id)  -  —  sin(«i)] 

K2  Kg  &2  K>2 

~~(f2  +  r)(l  +  T^){sin[2&2(d  -t)  +  kgt]  -  cos[2 k2(d  -t)+  kgt ]} 

2  &2  —  kz  k2  Kg 

-~(f2  +  f3)(l  -  ^){sin[2fc2(d  -t)-  kgt]  +  cos[2fc2(d  -t)- 
2  rC2  —  ^3  Kg 

+(cos  Kt+-  sin  Kt)fe  ~*~-f~)  sin[2A:2(d  - t)] 

AC  A^2  &3 

— cos  Kt  —  sin  Kt) ( f 2  +  )  cos[2A:2(d  —  t)] 

K^2  K>2  AC3 
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c  =  2(l-^)cos(kgt)-(l  +  ^)(l  +  ^)(^^)cos[kgt  +  2k2(d-t)} 
k2 


kgJ'~  ‘  k2'^k2-k3J 


k, 


k2  k2 


)  cos[fcst  —  2k2  ( d  —  £)] 


D  =  2^F  ~  ^  sin(M)  -  (!  +  ^K1  +  ^)(^  +  sin[V  +  2h(d  -  *)] 


)  sin[£;st  —  2&2  (d  —  t)] 


Fi  = 


)  sin  k2(d  —  t ) 


A: 

{[cos(fcst)  cos  k2(d  —  t)  —  -2-  sin (kgt)  f 
+i[7^  cos(&st)  sin  k2(d-t)  +  -£-  sin  (kgt)  cos  A^(d  —  t)]} 

=  Y\  Re  +  Im 

1.  k  8  ^ 

>2  =  - { [7^ (sin  kgt  —  —  cos  kgt)  +  —  (cos  Kt  —  —  sin  Kt)]  sin k2(d  —  t) 

2  k2  kg  k2  o  ' 

8  8 

+[(cos  nt  H —  sin  Kt)  —  (cos  kgt  +  —  sin  kgt )]  cos  k2(d  —  t)} 


Thus  the  field  solution  (23)  within  the  grating  may  be  expressed 

e,m  =  „{[(. 

+[(1  -  ~  \(l  ~  J")] e~ikgX  +  (COSKa;  +  ^smKx)}  (36) 


=  b(z)X(x ) 


where 


x(x)  =  [(i  +  -  I(I  +  i-)]e^ 


kg'C  +  iD  2' 


ikn 


(37) 


r/  fci  v  A.  +  %B  1 .  8  5  \ 

+[(1  ~  Tr)7TT7K  ~  o(2  ~  ~)]e  9  +  (cos kx  +  —  sin kx) 


ka'C  +  iD  2 


ikn 


K 


One  additional  relation  is  required  to  relate  the  decrease  in  energy  in  the 
wave  E0  given  by  (19b)  to  the  energy  radiated  from  the  guide  by  the  components 
Elx(z)  and  E3x(z).  Using  the  Poynting  vector  notation 


-LC6.^  +  5LCl.Cr  =  |(^e-2^) 


(38) 
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where  Aom  is  the  value  of  field  related  to  the  value  A  at  x  =  0  as  Aom  =  Ay(/c2  +  62)/2k 2 
Also,  rji  is  the  characteristic  impedance  of  the  medium  1  (771  oc  u/ni),  r]f  is  the  char¬ 
acteristic  impedance  of  the  waveguide  (r/f  oc  uj/n^  cos  6b)  and  dl  is  the  effective  width 
of  the  guide  from  which  the  power  in  the  input  wave  Eq  can  be  extracted.  This  re¬ 
lation  states  that  the  gradient  of  the  total  guide  power  in  z  is  equal  to  the  vertically- 
directed  power.  The  effective  width  of  the  guide  is 

d'  =  (d+\  +  -)  (39) 

0  7 

where  the  terms  d,  6 ,  7  are  due  to  the  mode  shape  which  comes  from  the  conventional 
calculation  of  power  inside  and  outside  of  the  guide  using  (7).  Equations  (34),  (35), 
(38),  (39)  are  solved  to  give 


2a  = 


kt  k2  B2(rlfA2  +  B2 
k2g(n2  -  k2)  k2  +  <52  d'  1 77  1  [C2  +  D2 

.Vf  rA2  +  B2  /t/ 2  .  \  1 


(40) 


where  the  definitions  of  7,  S,  k  are  given  by  (7).  Power  loss  is  defined  as  2a.  It  is 
evident  from  (38)  that  due  to  the  exponential  decay  of  input  power  with  z  then  2 a 
also  defines  the  diffraction  efficiency  as  2a  =  ay  =  (including  all  radiation  from 

•*171  9 

the  guide).  Actually aT  is  proportional  to  the  efficiency,  since  efficiency  as  a  unitless 
quantity  would  be  a?d!  but  we  will  refer  to  ax  in  general  as  the  efficiency. 


either  the 
as 


In  general,  one  would  be  interested  in  the  diffraction  efficiency  for  power  in 
upward  (Pu)  or  downward  direction  (Pd)-  These  are  expressed  respectively 


Oint 


Kq 


k2  sin2(^) 


k2(n2  —  k2)  k2  +  62 


ir2d! 


?!L(*±£.) 

7 ll{C2  +  D2) 


(41) 


which  is  and 


ad 


4  ±  M  ~  d l  *  (Y2  +Y2  )  (42) 

k2g(K2  -  k2g)  K2  +  62  7t2d!  773  C2  ~\~  D2  1Re  lW  1 

+Y2  +  Whz  +  BKsly, 

+  2  +  C2  +  D2  J 


which  is 


-§lL.  Thus  ax  =  otu  +  ad- 

*in 


3.2  Results  and  Discussion 

In  Fig.  4a  we  plot  the  dependence  of  a  upon  t,  the  grating  depth  for  both  the  TEq 
and  the  TE\  mode.  In  Fig.  4b  the  TE0  result  is  shown  for  two  different  waveguide 
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thicknesses.  The  wavelength  is  chosen  to  be  equal  to  the  grating  pitch  (i.e.  A/ n2  =  A) 
so  that  ideal  second  order  diffraction  is  obtained.  The  indices  and  other  parameters 
used  (indicated  in  the  figure)  are  chosen  to  coincide  with  the  calculations  of  Streifer 
and  some  of  his  data  points  are  also  shown  in  the  Fig  4a  for  the  TEq  calculation. 
(These  indices  are  typical  of  a  semiconductor  laser  waveguide  in  which  the  index 
differences  —  n\  and  n2  —  n2  are  small  for  a  weakly  guided  wave).  It  is  noted 
that  the  calculation  agrees  fairly  well  with  Streifer’s  although  his  is  slightly  smaller 
because  as  he  points  out,  he  has  not  included  the  contribution  of  the  E2  wave  in  his 
perturbation  analysis.  The  main  feature  of  the  result  is  the  resonance  that  occurs 
with  t.  At  t  —  0.14  yum,  constructive  interference  occurs  for  the  x  directed  wave  in  the 
grating  and  as  t  =  0.24  nm  the  interface  is  destructive.  Therefore  grating  thickness 
is  a  key  concern  when  optimizing  the  diffraction  efficiency. 

In  Fig.  5  the  variation  with  the  tooth  width,  w,  is  shown  and  as  expected 
a  maximum  is  found  for  a  slightly  asymmetric  tooth  shape  of  w/A  =  0.52.  As  shown 
this  behavior  also  agrees  fairly  well  with  Streifer’s.  Streifer  notes  that  the  maximum 
occurs  for  w/A  slightly  less  than  0.5  whereas  we  find  the  w/A  slightly  greater  than 
0.5.  Inspection  of  (40)  indicates  that  w/A  >  0.5  for  the  maximum. 

In  Fig.  6,  the  grating  thickness  is  fixed  and  the  waveguide  thickness  is 
varied.  It  is  noted  that  the  total  diffraction  efficiency  is  reduced  for  increasing  thick¬ 
ness  and  Streifer’s  calculation  is  also  indicated.  The  approximate  region  of  validity 
is  limited  to  a  ratio  of  t/d  =  0.7  as  discussed  later.  The  undulations  in  the  curve  are 
due  to  the  same  resonance  as  observed  in  Fig.  4. 

In  Fig.  7  the  index  72,3  is  varied  and  the  variation  of  diffracted  power  in 
the  upward  and  downward  direction  is  shown  together  with  the  total  radiated  power. 
Here  also  the  result  agrees  well  with  Streifer’s  exact  numerical  procedure. 

In  Fig.  8  the  variation  of  a  versus  t  shown  in  Fig.  4  cannot  be  valid  as  t 
approaches  d  as  mentioned  previously.  To  address  this  issue  it  is  useful  to  define  the 
effective  index  neff  for  the  propagation  of  the  wave  Eq  in  the  guide  in  the  presence 
of  the  grating.  We  use  the  relations 

0o  =  neffko  cos  0b  (43) 

0o  =  n^ko  cos  df,  (44) 

0o  —  nko  cos  9f)i  (45) 

to  ensure  a  constant  propagation  constant  in  the  all  regions  of  the  guide.  This  is 
illustrated  in  terms  of  a  ray  diagram  shown  in  Fig.  8a  which  shows  the  ray  AB  for 
the  guide  without  a  grating,  the  ray  AFC  for  a  composite  guide  consisting  of  both  n2 
and  n  regions  and  the  ray  AC  showing  the  final  ray  after  the  effects  of  propagation 
in  the  grating  are  included.  To  appreciate  the  effect  of  the  grating  in  this  diagram, 
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the  vector  (3q  is  indicated  in  the  vector  diagram  of  Fig.  8b.  The  grating  vector 
subtracts  from  /30  to  give  a  net  vector  of  length  /5X  according  to  (5)  as  shown.  From 
the  diagram  6b,  6b  1  can  be  found  from  geometrical  relations  in  order  to  obtain  the 
lateral  displacements.  Using  (43)-(45)  together  with  Fig.  8a,  the  effective  index  may 
be  found  as, 


neff  — 


ri2  cos  8b 
cos  6% 


n 2  cos  6b 

cos{tan~1  [  (d_t)  cotg+,cot  e61-]> 


n2  cos  6b 

( d—t )  cot  6b+t  cot  0b i 
y/d?+[(d—t)  cot  6b+t  cot  ]2 


(46) 


where  cot  db  = 


0o 


and  cot  = 


are  determined  from  Fig.  8a. 


In  Fig.  9  we  plot  the  variation  of  ne//  versus  the  grating  thickness  with 
different  ft  values.  The  value  of  neff  starts  from  at  t  =  0  and  slowly  increases  with 
increasing  t.  The  result  shows  the  perturbation  technique  is  still  reasonable  until  the 
grating  thickness  is  greater  than  about  0.7  of  the  waveguide  thickness  (±4%  errors). 
This  justifies  in  general  the  approximation  using  ne//  —  u2  as  discussed  earlier. 


3.3  Approximate  Second  Order 


Suppose  the  input  wavelength  does  not  coincide  with  the  grating  period,  so  that  the 
diffraction  is  either  somewhat  greater  or  somewhat  less  than  second  order  (non-zero 
values  of  6).  Then  the  differential  equation  is  given  by  (compare  with  (16)  and  (20)) 

+  (ky  -  0l)Ex  +  +  >4B  E0(z)  =  0  (47) 


where  now 


0i  = 


27 m2 
Ao 


cos  6b 


27 r 

T 


(48) 


which  now  includes  two  additional  terms  by  comparison  with  (20).  To  solve  this 
equation,  we  will  assume  that  E\  has  a  form  similar  to  the  ideal  second  order  result 
given  by  (36),  i.e.  we  assume 

El  =  b'(z)X'(x)  (49) 


where  U{z)  is  an  exponential  function  of  z  which  is  expected  to  be  similar  but  not 
necessarily  the  same  as  b(z)  as  given  by  (22).  Since  we  expect  the  same  decay  rate  in 
the  z  direction,  then  (47)  becomes 


cPEl 

dx2 


+  (kg)  Ei  +  k%  B  Eo(x,z)  =  0 


(50) 


where 

K  =  k'gr  +  ik’gi 


(51) 
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and  the  components  k'gr,  k'gi  are  found  by  splitting  k'g  =  \Jkln2  -  -  2ia/?i  into  real 

and  imaginary  parts  to  give 


k!  = 

Kgr 


kM—A  +i\/(^-/3f)2  +  4a^ 


1 1/2 

( 

/  2—2  o2\ 

-( 

k0n  -Pt) 

1/2 


(52) 


and 

*/ _ M (53) 

k$n*  -  (3?  +  (A^2  -  ft)  +  l°?0l\  *r 

since  a  -C  A,  fcon.  The  solution  procedure  now  follows  that  for  the  ideal  second  order 
and  equations  (27)-(32)  are  reproduced  with  the  modification  that  kg  is  replaced 
wherever  it  appears  by  k'g  and  that  U(z)  now  has  real  and  imaginary  parts  given  by 


b'(z)  =  br  —  ibi  (54) 

k2n2(K2-k2n2  +  p2)  ,  2 a/h  A  c-~ 

~  (k2  —  koH2  +  0i)2  +  4a2/?2  («2-/c§n2  +  /?2)2  +  4a2^J  * 


The  results  are  similar  to  (34)  and  (35)  with 


r  *h'A' +  iB'  r 
Cl  ~ 2b  aTi&’  Ct 


(55) 


It  will  be  appreciated  that  the  solution  of  the  linear  equations  (22)-(32)  with  these 
changes  is  tedious  by  hand  but  that  a  simple  matrix  solution  is  possible  and  readily 
solved  by  conventional  programs  to  obtain  relatively  quick  results.  On  the  other  hand, 
closed  form  solutions  for  A',  B',  C',  D\  b',  Y[  and  can  be  found  straightforwardly 
and  then  results  can  be  plotted  in  about  3  seconds  on  a  PC  platform.  Some  of  these 
results  are  given  in  the  appendix  C. 


In  order  to  express  the  diffraction  efficiency  in  a  similar  way  to  the  ideal 
second  order  case,  it  is  assumed  that  a  can  be  used  for  a0  in  (52)-(53).  This  expedient 
is  justified  as  we  show  later.  Then  the  expressions  for  au  and  ad  as  given  by  (41)  and 
(42)  can  be  used  with  A,  B,  C,  D,  Yj  and  Y2  replaced  by  their  primed  counterparts. 

These  results  are  illustrated  in  Fig.  10  which  shows  the  variation  in  the 
diffraction  efficiencies  a  (total),  au  (upwards)  and  ad  (downwards)  as  a  function  of  A. 
Also  shown  on  this  plot  using  the  points  is  a  comparison  of  a  calculated  by  using  c*o 
for  a  in  (52)-(53)  versus  a  totally  self-consistent  calculation  and  there  is  no  difference 
in  these  results.  The  value  of  a  is  seen  to  go  to  a  minimum  at  A  =  2n//3o  which 
is  the  ideal  second  order  condition  and  to  rise  for  a  smaller  or  larger  values.  This 
implies  power  escaping  from  the  guide  and  propagating  in  both  x  and  z  directions 
with  the  z  component  increasing  as  we  move  away  from  the  second  order  condition. 
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This  is  also  explained  by  the  presence  of  k'gi  which  implies  an  x  dependent  decay  in 
the  wave  CxeiklX.  The  decay  in  energy  is  accounted  for  by  the  z  travelling  wave  with 
propagation  constant  (3 i. 

Now  for  the  case  that  the  diffraction  is  not  perfectly  second  order,  there  is 
a  component  of  diffracted  power  flow  both  in  the  x  direction  and  in  the  z  direction. 
If  A  <  A,  the  m  =  1  component  of  the  field  has  z  diffracted  power  primarily  in  the 
same  direction  as  the  incident  wave  and  if  A  >  A,  the  component  is  primarily  in 
the  opposite  direction  to  the  incident  wave.  Therefore  the  direction  of  the  m  —  1 
diffracted  wave  makes  an  angle  with  the  vertical  direction  of  6X  where 

tan  6\  =  =  —  cos  - y  (56) 

k\  rii  n  iA 

Consequently  the  x  axis  in  Fig.  10  also  represents  the  angular  deviations  of  the  total 
emitted  energy  with  respect  to  the  vertical  direction  (ideal  second  order)  and  this  is 
indicated.  The  discussion  here  is  limited  to  values  of  9\  which  are  actually  limited 
by  (24)  or  (26)  since  h  and  k3  must  be  real  in  order  to  propagate  power  out  of  the 
waveguide.  Therefore  using  (24),  (26)  and  (56)  the  limitations  on  9\  are  given  by 

9lu  <  tan'1  ^  ,  9ld  <  tan"1  -  (57) 

Til  Til 

These  limiting  values  are  about  ±45°  using  the  parameters  in  Fig.  10a. 

Since  these  values  depend  on  the  index  difference  in  the  grating  then  for 
weak  gratings  it  may  be  the  case  that  9\u  >  63rd  (see  Fig.  2b).  For  these  situations 
there  is  no  distinct  regime  of  grating  pitch  without  diffraction.  Rather  there  is  a  con¬ 
tinuous  coexistence  and  overlap  of  the  approximate  second  and  third  order  diffraction 
components. 


4  Analogy  to  coupled  mode  theory 


Coupled  mode  equations  are  used  to  determine  the  guided  wave  solutions  to  the  first 
order  waveguide  diffraction  which  would  be  the  E0  and  E2  solutions  in  this  paper.  The 
coupled  mode  approach  is  to  represent  the  grating  diffracted  waves  as  conventional  T E 
or  TM  waves  and  then  substitute  into  the  perturbed  Maxwell  equation  to  determine 
the  wave  envelope  shape  (i.e.  the  growth  or  decay).  The  present  approach  to  second 
order  diffraction  is  analogous  to  the  coupled  mode  approach  in  that  the  first  order 
waves  are  represented  in  the  same  way  in  order  to  predict  the  radiating  wave. 
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5  Summary 


The  diffraction  efficiency  of  a  second  order  grating  has  been  calculated  by  considering 
first  the  ideal  second  order  situation  of  emission  at  90°  to  the  guide  and  then  arbitrary 
angles  of  emission  as  separate  cases.  A  closed  form  result  has  been  found  based  on 
the  approximation  of  a  single  value  of  propagation  constant  (/?o)  for  the  guided  wave. 
The  approach  maintains  accuracy  by  including  both  forward  and  reverse  first  order 
diffracted  waves  as  perturbation  terms.  The  magnitude  of  the  composite  wave  is  not 
important  since  it  cancels  out  in  the  determination  of  a. The  calculations  are  compared 
with  the  numerical  results  of  earlier  work  and  it  is  shown  that  the  analytical  or  quasi- 
numerical  simplified  approach  is  a  useful  tool  for  iterative  device  design.  An  important 
aspect  of  this  work  is  that  the  solutions  are  valid  for  strong  or  weak  gratings  because 
no  approximations  have  been  made  which  depend  on  a  small  index  change. 
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Appendix  A 

Particular  Solution  of  Wave  Equation 

The  x  dependent  differential  wave  equation  in  the  grating  is  obtained  from 
(19b)  and  (20)  as 

^  +  ^,  =  -*{8^-*=  (Al) 

The  complementary  solution  (RHS=0)  is 

Egc(x)  =  C2eik°x  +  Cze~ik*x  (A2) 


the  particular  solution  is 
5  A*  f 


Enn  =  ki  B  e1' 


J9P  ~  0  “  ~  I  -2 ikQ 


7 


A0(x)e~azeik9X 


■dx 


k2  ~  _  6 

=  (  o  °  ,o)  B  A-e  a~(cosKX  +  -  sin kx) 

Kr  —  kg  K 


(A3) 


where  —2 ikg  is  the  Wronskian  and  A0(x)  in  (7)  has  been  used  for  the  grating  region 
to  represent  TE  modes.  Then  the  total  solution  is 


Egc(x )  =  C2eik*x  +  C3e~tK°x  + 


ik„x  ,  k0 


.  , .  B  Aye  az (cos  kx  H —  sin  kx')  (A4) 

K“  —  kg  "  K 
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Appendix  B 

Sinusoidal  Grating 

For  the  sinusoidal  grating,  the  coefficients  in  (1)  are  expressed 

— !!ll<S2n[gcos-1(^  —  1)]  0  <  x  <  t 

Bq  =  <  nq 

0  r  0,  x  t,  q  0 


(Bl) 


Then  equation  (13)  is  modified  to 

-  0}Ei  +  k\f?Ex  -  kl— — —  sin(q cos-1  ( ^  -  l))£o(x)  =  0  (B2) 

oxz  nq  t 

Then  the  solution  for  Eg  in  the  grating  is  written  as 

Eg(x,  z)  =  C2eik°x  +  C3e~ik°x  +  Egp  (B3) 

and  the  particular  solution  can  be  obtained  by  following  the  procedure  in  Appendix 

A, 

E„  =  b  =  H  (nl  -  sinfe(x  -  |)j  A0ge^  (B4) 

where  Jx  is  the  Bessel  function  of  the  first  kind. 

The  remaining  procedure  and  results  are  similar  to  those  of  the  rectangular 

grating. 
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Appendix  C 


Solution  of  the  linear  equations  (22)- (32)  for  the  approximate  second  order  diffraction 
condition  yield  the  following  results  for  the  parameter  A',B',C',D'. 

A  =  {■p——r)  cos2k2(d  —  t)[(l  —  8)  cosKt  +  (—  +  k)  sin  Act] 
k2  —  «3  K 

—  [(1  +  8)  cos /ct  +  ( - k)  sin  Act] 

AC 

1  8k"  8  k' 

+5[(1  -  cos^  +  “pHW 

*[1  +  T--(z—r)(1  -  T-) cos  2*a(d  - 1)  +  (f* +  sin2A ^(d  -  <)] 

A?2  rC2  —  /C3  rC2  A?2  ““  ^3  & 2 

«...  « 


+2[(1  “  kfTkf'>^t ~  cos*»<1  exp(_<:’4) 

&2  +  ^3^  _  ^ 

k2 

8k' 


k' 

p 

l/c2 


♦It*-  +  ( 


,&2  +  ^3  n  kg 
■  ^3  h 


)(1  -  ^)sin2 k2(d-t)  +  (y  '  7/ ) ^  cos2k2(d  -  t)] 


k2  -  k3 

1  XZ-"  A/-' 

+  2[(1  +  C°S/C^  +  feg>+V  Sin^]  eXp(^t} 
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Fig  1  A  slab  waveguide  with  a  grating  of  thickness  t 
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Fig.  2a  Bragg  diffraction  scheme 
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Fig.  3  Schematic  of  the  incident  wave  producing  transmitted  wave, 
reflected  wave  and  diffracted  wave  due  to  the  grating 
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Fig.  4a  Diffraction  efficiency  versus  grating  thickness  for  a  rectangular  grating. 
Discrete  points  are  data  of  Streifer[10].  •  TE0  mode  ■  TE,  mode 
Note  that  the  individual  points  (from  Streifer's  paper)  are  not  for  the 
ideal  second  order  case.  Therefore,  even  though  these  should  be 
larger  than  our  ideal  second  order  calculation,  it  is  not  because 
Streifer  omits  E2  from  his  perturbation  (see  text). 
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GRATING  THICKNESS,  t  (jxm) 

Fig.  4b  Diffraction  efficiency  versus  grating  thickness  for  a  rectangular  grating 
with  different  waveguide  thicknesses. 


TOOTH  WIDTH,  w/A 

Fig.  5.  Diffraction  efficiency  versus  tooth  width  for  a  rectangular  grating 
Discrete  points  are  calculations  of  Streifer[10]. 


DIFFRRACTION  EFFICIENCY,  a  (cm  ') 


WAVEGUIDE  THICKNESS,  d  (fim) 

Fig.  6  Diffraction  efficiency  versus  waveguide  thickness  for  a  rectangular  grating 
Discrete  points  are  calculations  of  Streifer[10]  for  the  TE0mode. 


DIFFRACTION 


Fig.  7  TEq  mode  diffraction  efficiency  versus  n3  for  a  rectangular  grating. 
Discrete  points  are  calculations  of  Streifer[10]. 


Fig.  8a,  8b  Schematic  of  total  internal  reflection  for  a  composite  waveguide 
and  wave  vector. 
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9  Normalized  effective  index  (with  respect  to  n2)  versus 
grating  thickness  with  Pi  as  parameter. 


DIFFRACTION  EFFICIENCY,  a  (cm  ') 


Fig.  10a  Diffraction  efficiency  versus  grating  pitch  for  approximate  second-order 
rectangular  grating.  Discrete  points  are  solutions  obtained  with  ao  to 
represent  a  in  /  PidEi/cfz  (=/  PiOtEi)  in  (47). 
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DIFFRACTION  EFFICIENCY,  a  (cm  ') 


Fig.  1 0b  Diffraction  efficiency  in  the  upward  and  downward  directions  versus 
grating  pitch  for  the  approximate  second-order  rectangular  grating. 
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Appendix  B.  Diffraction  Into  Corrugated  Waveguide  From  Normally  Incident 
Radiation 


Diffraction  Into  a  Corrugated  Waveguide  From 
Normally  Incident  Radiation 


G.  W.  Taylor  and  C.  Kwan 
Electrical  and  Systems  Engineering 
University  of  Connecticut 
Storrs,  CT  06269 


Abstract 

The  diffraction  of  a  plane  wave  normally  incident  upon  a  corrugated  waveg¬ 
uide  with  a  second  order  grating  pitch  is  examined.  Using  a  rectangular  grating 
the  efficiency  of  diffraction  into  the  first  order  components  in  the  form  of  guided 
waves  is  determined.  Using  simple  approximations,  closed  form  expressions  are 
found  from  which  the  dependence  of  efficiency  upon  grating  thickness,  waveg¬ 
uide  thickness,  grating  pitch  and  indices  are  predicted.  The  approach  provides 
a  useful  simulation  tool  to  predict  the  efficiency  of  coupling  into  a  waveguide 
from  an  external  or  internal  orthogonal  source. 
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I  Introduction 


Coupling  into  or  out  of  a  waveguide  by  the  use  of  a  grating  has  been  studied  by  many 
researchers  [1]  [2]  [3]  [4]  as  a  simple  method  to  transfer  free  space  signals  to  waveguides 
in  optical  integrated  circuits.  In  most  cases  [2]  [4]  the  input  light  is  propagated  at 
an  angle  to  the  guide  to  optimize  the  coupling  since  the  tangential  input  component 
is  launched  more  easily  into  the  guide.  However  propagation  normal  to  the  guide 
is  also  important  as  a  means  to  couple  from  a  vertically  emitting  structure  into  a 
waveguide.  With  the  emergence  of  the  vertical  cavity  laser  as  a  preferred  embodi¬ 
ment  to  an  integrated  light  source,  the  use  of  a  normally  incident  diffraction  grating 
becomes  an  interesting  mechanism  to  create  an  optical  source  with  high  right  angle 
mode  output  power[5][?].  Another  new  important  application  for  diffraction  from 
normally  incident  light  to  laterally  propagating  light  is  the  detection  of  infrared  light 
using  intersubband  transition  [6].  It  is  therefore  of  considerable  interest  to  develop 
predictions  of  the  efficiency  of  the  diffraction  process  in  terms  of  the  waveguide  and 
grating  parameters.  The  existing  mathematical  descriptions  [1][3][4][7],  which  in  gen¬ 
eral  describe  off  normal  incidence  situations,  are  numerical,  while  our  objective  here 
is  to  obtain  an  analytical  prediction. 

Previously  [8]  we  have  reported  an  approach  for  describing  diffraction  from 
the  guide  into  the  superstrate  or  the  substrate  which  enables  a  closed  form  solution. 
The  same  techniques  will  be  applied  here  to  characterize  the  diffraction  in  the  reverse 
direction.  For  this  derivation  a  symmetrical  grating  will  be  assumed. 
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II  Diffraction  into  guide 


In  the  previous  paper  [8]  Streifer’s  method  [1]  was  used  to  describe  the  problem 
of  ideal  second  order  diffraction  of  a  guided  wave  into  a  wave  radiating  at  right 
angles  to  the  guide.  Here  we  adapt  the  same  approach  to  the  diffraction  of  a  wave 
normally  incident  upon  a  corrugated  guide,  into  a  guided  mode  as  illustrated  in 
Fig.  1.  The  slab  waveguide  is  corrugated  to  a  thickness  t,  an  arbitrary  point  in  the 
corrugation  is  chosen  for  the  origin  of  the  z  axis,  the  x  axis  origin  is  chosen  at  the  top 
of  the  corrugation  and  x  increases  positively  downwards.  As  previously,  the  periodic 
variation  of  the  index  within  the  grating  is  represented  by  the  Fourier  series 


n](x,z)  =  r?(x)+  Y.  Bt(zy2,,'/A  (1) 

q——oo  &:  t^O 

where  fi  is  the  average  index  over  one  grating  period,  Bq  are  the  Fourier  coefficients, 
and  A  is  the  grating  period  as  illustrated  in  Fig.  1.  For  the  rectangular  grating,  which 
is  our  main  interest  here,  these  coefficients  have  been  evaluated  as  (see  Streifer  [1]) 


Bq 


(r4-nl)sin{^W) 


7 rq 


0 


0  <  x  <  t 

x  <  0,  t  <  x,  q  0 


(2) 


where  w  is  the  grating  tooth  width  (see  Fig.  1).  By  restricting  the  discussion  to 
include  only  the  first  and  second  diffraction  orders,  the  index  in  the  grating  region 
corresponding  to  (1)  is 


n2g(x,  z,  q)  =  n2  — 


{n: 


2  —  n?} 


i) 


f7rqw 


qix 


sin(—^~)e  A  q  =  ±1 


(3) 


where  the  average  n  is  given  by 


n2  =  {wn\  +  (A  —  w)n\}/A  ,  0  <  x  <  t 


(4) 
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Overall,  the  index  variation  is  summarized  by  region  as 


n2(x,z) 


n\  x  <  0 

n2(x,  z)  0  <  x  <  t 

n2  t  <  x  <  d 

n\  d  <  x 


(5) 


III  Second  order  grating  calculation 


Because  we  are  restricting  our  discussion  to  the  case  of  second  order  diffraction,  then 
it  was  shown  [8]  [9]  that  the  problem  is  simplified  because  we  need  to  consider  only 
the  first  and  second  order  diffracted  waves.  In  Fig.  1  the  normally  incident  radiation 
( 9  =  0)  corresponds  to  ideal  second  order  diffraction  and  other  waves  for  9  ^  0, 
represent  approximate  second  order  diffraction.  By  representing  the  total  E  field  as 
the  sum  of  the  diffracted  waves  in  the  forward,  the  reverse  and  the  surface  normal 
directions  and  substituting  into  the  two  dimensional  wave  equation,  it  was  shown 
that  the  diffraction  is  represented  by  three  differential  equations  (D.E.’s)  which  are 

-  $E0  +  k2n2E0  +  2i(30^  +  k2  (B1  Ex  +  B2E2)  =  0,  m  =  0  (6) 

ox1  oz 

-  0iEx  +  k2n2E1  +  +  k2  (B_!  E0  +  Bx  E2)  =  0,  m  =  1  (7) 

^  -  P2E2  +  k2n2E2  +  +  k2  (B. i  E1  +  B.2  E0)  =  0,m  =  2  (8) 

It  is  noted  that  the  coefficients  q  =  —  1  and  +1  in  (2)  determine  the  coefficients  for 
the  waves  with  m  =  0  and  2  (i.e.  waves  Eq  and  E2)  respectively.  These  will  be 
represented  by  B- 1  =  Bi  =  B  = 

These  three  diffracted  waves  are  indicated  in  Fig.  1  and  the  variation  of 
the  Eq  and  E2  components  in  the  z  direction  is  indicated.  These  components  are  the 
first  order  diffracted  waves  and  are  growing  in  the  +z  and  —  z  directions  respectively 
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due  to  diffraction  into  the  guide  from  the  normally  incident  wave.  It  is  interesting  to 
note  however  (and  specified  later)  that  second  order  diffraction  is  required  to  produce 
Eq  and  E2,  i.e.  they  would  not  exist  with  surface  normal  radiation  for  a  first  order 
diffraction  grating. 

The  generalised  Bragg  condition  which  relates  the  propagation  constants 
of  the  various  diffraction  orders  m  is  /3m  =  /30  -  For  the  ideal  second  order  case 
of  /3i  —  0,  it  was  shown  that 

00  =  -02  (9) 

and  since  m  =  1  for  this  case,  then  we  have  0q  =  jL.  Then  the  E\  D.E.  (c.f.(7)) 
becomes 

pfirp  _ 

+  k^E,  +  k20  B(  E0  +  E2)=  0  (10) 

and  we  also  note  that  for  the  symmetrical  grating,  we  have  E0  =  E2  (the  case  of 
a  .symmetry  will  be  considered  elsewhere).  The  Ex  wave  exists  in  the  grating  region  and 
results  from  the  wave  Coe~iklX  incident  on  the  waveguide.  Due  to  multiple  reflections 
at  the  interfaces  0,  t  and  d,  there  will  be  several  other  normally  propagating  waves 
within  the  structure  which  are  shown  in  Fig.  1. 


In  regard  to  the  waves  Eq  and  E2,  we  utilize  the  model  of  the  previous 
paper  for  diffraction  out  of  the  guide.  Thus  E0  and  E2  are  assumed  to  be  guided 
waves  traveling  in  the  ±z  directions  in  the  guide  and  represented  by 

E0  =  Ao{x)eaze-^z ,  E2  =  Ao{x)e~az eipoz  (11) 


where  fio  is  the  propagation  constant  of  the  waveguide  and  a  is  the  growth  constant 
of  the  wave.  These  waves  are  assumed  to  consist  of  TE  modes  in  the  guide  with 
Aq(x)  determined  by 


Tlo(rr)  =  < 


Aye6x 

Ay  cos  kx  +  By  sin  kx 

(Ay  cos  nd  +  By  sin  Kd)e~l('x~d) 


x  <  0 
0  <  x  <  d 
d  <  x 


(12a) 
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where 


Ay  is  the  value  of  Aq(0), 


k  =  n2ko  sin  9b, 

8  =  (n2  —  n\)  /cq  —  n\kl  sin2  9b 
7  =  [(n2  nt)  kl  ~  n2k0  sin2  Qb 


1/2 

1/2 


(12b) 


and  9b  is  the  bounce  angle.  As  mentioned  earlier  these  waves  are  assumed  to  grow 
in  the  ±z  directions  due  to  the  diffracted  energy  from  the  input  wave  and  a,  the 
rate  of  growth  will  be  determined  from  the  solutions.  The  presence  of  the  Eq  and  E2 
waves  represent  an  energy  loss  mechanism  for  the  wave  E\  in  the  grating  region  and 
(11)  provides  a  representation  for  the  perturbation  terms  in  (10).  Solving  (10)  using 
(11)  with  Eo,  f?2  =  0  to  find  the  complementary  solution,  using  (12a,  b)  to  find  a 
particular  solution  and  then  combining  the  two,  a  general  solution  (see  Appendix  A) 
is  found  to  be 


Egx(x,  z)  =  C2etkgX  +  Cze  lkgX  +  6(cos  kx  +  —  sin  kx) 

K 


(13) 


where 

01.2  _ 

B  Av*+aZ  (14) 

This  solution  introduces  a  factor  of  2  in  (10)  because  Eq  =  E2  for  the  symmetrical 
situation.  It  is  noted  in  this  case  that  the  presence  of  Eq  and  E2  in  (10)  represents  a 
small  loss  term  in  the  general  Ei  equation  rather  than  a  driving  force  in  the  out  of 
the  guide  case,  and  it  is  therefore  expected  that  Eq,  E2  E\.  Thus  E\  exists  even  if 
Eo,  E2  =  0,  But  for  the  Eq  and  E2  equations  (6)  and  (7),  the  perturbation  term  Ex 
is  a  driving  force  and  the  waves  Eq  and  E2  would  not  exist  without  it.  It  noted  that 
both  Ay  and  a  have  been  introduced  as  unknowns  which  depend  on  the  input  wave 
Co¬ 


lt  is  noted  that  Egx(x,z)  represents  the  E\  wave  solution  in  the  grating 
region.  In  the  work  that  follows  we  continue  to  use  Egc  for  the  x  directed  wave  in 
the  grating  and  use  E\x,  E2x,  and  Egx  for  the  x  directed  waves  in  the  other  parts  of 
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the  guide  which  represent  collectively  the  E1  wave.  They  are  described  by  region  as 
shown  in  Fig.  1  and  in  total  the  waves  may  then  be  written 


Elx(x,  z)  =  C0e~iklX  +  CieiklX  ,  kx  =  ( k20n\  -  pf)1/2  =  k0ni  (15) 

Egx(x,  z)  =  C2eik*x  +  C3e~ik°x  +  Egp  ,  kg  =  ( k2n 2  -  p2)1/2  =  k0n  (16) 

E2x(x,  z)  =  C4eik2^  +  C5e-ik2(d~x>  ,  k2  =  ( k2n\  -  02)1'2  =  k0n2  (17) 

E3x(x,  z)  =  C6e~iks(x~d)  ,  k3  =  (kin2  -  $)1/2  =  k0n3  (18) 


in  which  Egp  is  the  grating  dependent  contribution  as  described  by  (13)  and  fa  —  0 
as  given  by  (9). 

We  then  perform  the  matching  of  these  equations  in  both  magnitude  and 
derivative  (continuity  of  Hz)  at  the  interfaces  x  =  0,  x  —  t  and  x  =  d.  In  line  with 
the  previous  derivations,  the  matching  produces  the  equations 

Co  +  C\  =  C2  +  C3  +  b  ,  x  =  0  (19) 

— ik\Co  ik\C\  —  ikgC2  —  ikgC3  +  6  ■  b  ,  x  —  0  (20) 

C2eik^  +  Cze~ik9t  +  6 (cos  Kt+-  sin  Kt)  =  C4eifc2(d'°  +  C^2^  ,  x  =  t  (21) 

tz 

ikg(C2etk,9t  —  C3e~lk9t)  —  £>(/csin/ct  —  6  cos  Kt)  =  — —  C3e~tk2^d~^)  ,  x  =  t 

(22) 

C4  +  C5  =  C6  ,  x  =  d  (23) 

— ik2C4  +  ik2C$  =  — ik3Ce  ,  x  =  d  (24) 

When  these  equations  are  solved  we  obtain  the  following  results  for  Ci,  C2,  C3,  Cq, 

_  +  ^  +  (!-  ife)Qi  +  (1  +  + 

1  C1  +  ^)°i  +  C1  “  ^)°2  °  C1  +  *J)ai  +  (!  ~  §)a2 

=  qCo  +  r2b, 

c2  =  ho  -  h + (1 + y  n + ki + y  -  a + j-)]s  (26) 

Z  Kg  Kg  £  Kg  L!\*g 

=  t3Cq  +  r4b 
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c3  =  2^1  +  F^  +  ^1_F^C,°+2^1_  ’k^7'2~(1~ik~^b  ^ 

Zi  fvg  t Vg  Zj  rog  6/vp 

=  r5C0  +  r6b 


C6  =  (r3  •  g^-Md-O+M  +  r&  .  e*[fc2(d-0+M)(70  _|_ 

(Ii(l  _  hl')e*[-k2(.d-t)+kgt]  _|_  ^6^  _|_  ^£^ei[fc2(4-«)+A:st]  _|_ 

2  &2  2  k2 


(28) 


5  1 

(cos  Kt  +  -  sin Kt)  cosh[A:2(d  —  t)]  —  —  (k  sin  Kt  —  6 cos  Kt)  sinh^d  —  t)]}b 

K 

=  r7C0  +  r&b 


where  the  following  functions  have  been  defined 


0l  =  [(1  +  ^ )x  -  (1  -  tl)e-^2(d-t)ykgt  (29) 

a2  =  [(1  -  !%)x  -  (1  +  tl)e-^2(d-t)^e-ikgt  (30) 

K>2 

c  o  • 

a3  =  — 2(cos nt-\ — sin/ct)(^+e_2lfc2^-^)+— (/csin/ct— <5cos/ct)(— x+e-2lfc2^-^)  (31) 
K  k>2 


l  +  fc3/fc2 

1  -  h/k2 


(32) 


In  order  to  determine  a  we  apply  the  conservation  of  power,  which  de¬ 
mands  that  the  difference  between  power  input  normal  to  the  guide  and  the  power 
output  normal  to  the  guide  should  equal  the  gradient  in  z  of  the  guided  wave  power 
which  is  stated 

2T^,ALe+2°’d,) = ^Co '  c°'  -  kCi  ■  c‘  -  kc°  ■ Ci  (33) 

Here  2  in  the  left  hand  side  of  (33)  indicates  the  outputs  in  both  +  x  and  -  x  directions 
and  Aom  is  the  maximum  value  of  field  of  the  guided  wave,  given  by 

(34) 

where  d!  is  the  effective  guide  thickness  d!  =  d  +  |  ^  and  Ay  was  defined  by  (12a), 

(12b).  As  in  the  previous  work  [8]  we  assume  that  the  index  of  the  propagating  wave 
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in  the  corrugated  waveguide  is  n 2  which  was  shown  to  be  a  good  approximation  until 
the  grating  thickness  t  becomes  a  substantial  fraction  of  the  guide  thickness  d. 

To  proceed  further  with  the  solution  we  need  to  express  Aom  (or  equiv¬ 
alently  b  as  per  (14))  as  terms  proportional  to  Co,  the  input.  Then  all  terms  in  b 
on  both  sides  of  (33),  (i.e.  in  Aom,  C\  and  C&  )  can  be  cancelled  out.  This  can  be 
accomplished  by  solving  (6),  the  D.E.  for  Eq  corresponding  to  (11)  and  using  the  so¬ 
lution  (13)  with  C2  and  C3  given  by  (26)  and  (27)  to  represent  Ex  in  the  region  of  the 
grating.  After  matching  this  solution  at  all  boundaries,  the  relation  £o(0)  =  Ay  can 
be  used  to  determine  b  in  terms  of  Co. This  procedure  will  be  implemented  shortly. 


However,  in  order  to  determine  a  alone  this  procedure  is  not  initially 
required  if  we  note  that  in  the  absence  of  the  grating,  we  would  have  b  =  0  and  then 
(33)  reduces  to 

111 

(35) 


111 

— c0-c0  =  — c1-Ci  +  — a6-c6 


It  is  therefore  evident  by  inspection  of  (33)  and  (35),  that  only  the  terms  in  b  on  the 
right  side  of  (25),  (28)  are  required  and  so  using  (14),  (25)  and  (28)  in  (33)  we  obtain 


\_,k2  +  <52w«2  -  k 2 
V2 


=  ZTr 2  •  r2  +  f  rs  •  r< 


1_ 

V3 


(36) 


2k2  /y  k20B  '  Vi 
where  r2  and  r8  are  defined  by  (25)  and  (28).  It  is  noted  in  performing  the  complex 
conjugations  to  obtain  (36)  from  (35)  and  (33),  the  terms  in  C0  and  b  in  (25)  and 
(28)  represent  waves  travelling  in  orthogonal  directions  and  therefore  cross  products 
are  not  required.  The  solution  for  the  wave  growth  constant  a  is  obtained  from  (36) 
as 

01  =  (cmzKz-o — n)  i\-r 2  ■  r2  +  -rs  •  (37) 


Performing  the  complex  conjugations  a  positive  real  value  for  a  is  obtained  and  is 
expressed  in  terms  of  the  constants  A  —>  H  where 


A  +  iB 

~cTw 


=  r2, 


E  +  iF 

gTTh 


=  r8 


(38) 
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The  expressions  for  A,  B,  C,  D,  E,  F,  G,  H  are  determined  in  the  Appendix  B. 


Comparison  of  the  solution  of  A  — >  H  with  the  solutions  for  diffraction 
out  of  the  guide  [8]  reveals  that  the  result  obtained  for  a  is  identical  for  the  two  cases. 
This  may  also  be  seen  by  inspection  of  (19)  (24)  since,  if  Co  is  reduced  to  0,  this  set 
becomes  identical  to  that  for  diffraction  out  of  the  guide.  In  some  sense,  this  result 
is  expected  as  the  natural  reversibility  of  light  propagating  through  a  passive  system 
since  the  decay  in  one  direction  becomes  the  growth  in  the  opposite  direction  of 
propagation.  However,  this  reversibility  does  not  apply  to  the  efficiency  as  described 
below. 


In  the  previous  derivation  of  diffraction  out  of  the  guide,  the  decay  con¬ 
stant  of  the  input  wave  times  the  effective  guide  width  was  the  diffraction  efficiency. 
For  diffraction  into  the  guide,  this  is  no  longer  true  due  to  the  significant  transmission 
through  the  guide.  Furthermore,  since  the  input  normal  to  the  guide  is  constant  with 
position  and  yet  the  wave  diffracted  into  the  guide  is  growing  with  distance  then  the 
overall  efficiency  would  appear  to  increase  with  increasing  z.  This  is  erroneous  of 
course  because  the  mode  growing  in  the  guide  would  also  diffract  back  out  of  the 
guide  to  a  radiating  wave.  Therefore  to  determine  the  efficiency  of  diffraction  into 
the  guide  as  an  independent  mechanism  on  a  per  unit  length  basis  at  an  arbitrary 
position  one  should  consider  z  =  0,  since  by  definition  this  result  could  not  be  affected 
by  any  diffracted  waves  for  z  <  0. 


V  = 


dPrw 

SPin 


(39) 


where  8Ptw  is  the  TW  power  generated  in  an  increment  dz  and  8Pin  is  the  power 
incident  in  the  corrugation  in  an  increment  dz.  From  the  power  components  in  (33), 
this  is  determined  as 


r]  =  2  ad' 


Vi 


A  A* 


r)2  CqCq 

A  expression  similar  to  (40)  can  also  been  found  in  another  publication[10]. 


(40) 
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An  issue  of  importance  is  how  to  use  the  two  independent  results  of  diffrac¬ 
tion  into  the  guide  and  diffraction  out  of  the  guide  in  a  practical  application  such  as 
determining  the  length  of  guide.  In  general  one  must  write  a  differential  equation  in¬ 
corporating  both  effects  separately  in  addition  to  other  losses  or  sources  in  the  guide. 
In  the  absence  of  these,  the  two  results  can  be  readily  combined  to  find  a  length 
dependent  efficiency  r](L). 


To  determine  rj  we  require  Ay  and  this  means  a  solution  of  (6)  is  required 
in  all  regions  of  the  guide.  We  therefore  write  the  equations  for  the  top  cladding  and 
the  grating  as 


Eta  =  ,  k[  =  -  ;3;)1/2 


(41) 


and 


*(.,  *>  -  A*"  +  (42) 

Kg  Kg  Kg  Kg 


Noting  that  the  constants  can  be  rewritten  as 


0o  =  fcon2 cos2  fa,  k'g  =  kin2  -  01 4-  2 ia/30,  k2g  -  k'2  =  01 

and  if  we  ignore  the  term  2ia/3o  since  a  Po  then  (42)  becomes 

E0g(x ,  z)  =  D2eik>  +  D3e-ik*x  +  B\C2eik*x  +  C3e~ik*x)  (43) 

where  B'  =  2  ^~rr-  For  the  waveguide  core  and  lower  cladding,  the  equations  are 

77.2  COS  pj 

Em(x, z)  =  DieLtVd->  +  ,  k'2  =  (*02n|  -  ;52l,/2  (44) 

Eta(x,  z)  =  Dee-^-^  ,  It ,  =  (fcgnf  -  /302)1/2  (45) 


Equations  (41)-(45)  are  matched  at  each  interface  in  both  magnitude  and 
derivation  as  implemented  for  the  Eiwave.  Then  it  is  possible  to  express  Dx  in  terms 
of  C2  and  C3 


£>i  = 


yi  •  C2  +  ^2  •  Cz 
yo 


(46) 
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where 


yo  =  I(i  +  §)(i  +  h.)x>eiw2(.d-t)+ik>gt  +  I(i  _  M)(i  _  b.y^2{d-t)-ik'gt  ^ 

2  kg  ^  2  kg 

+ ^)(i  -  !)e^  -  ^(i  -  ^)(i + l)e_^ 

^  A/j  A/2  ^  A/^  /L'2 


^ 9\J,kQt  ( i  |  ^9\„ikgt  t  i2k'2(d-t) 


yi  =  B'{(  l-^K^-a  +  TfJc^Vc 

+i(l  +  #)e^*[(l  +  %)x'eW-Q  -  (1  -  §)] 

Z  K  K  2  K'2 


+J(1  -  |)e-^[(l  -  §)^(4W)  -  (1  +  #)]} 


fc. 


A:: 


& 


(48) 


y2  =  £'{(l  +  jj)e~ikst  -  (1  -  ^.)e~ik^x' ei2k'^d~t) 
k<2  k2 

+|(1- j*)<*‘[(l  +  -(!-§)] 

Z  /v,,  (V9  A/9 


k'2 

K 
k'2 


+J(1  +  |)^[(1  -  %)x'eik>W  -  (1  +  &)]} 

/Co 


6'  = 


2k$B 

(*2  ~  k?) 


,  x  = 


i  +  ^3/^2 
1  -  **/*£ 


(49) 


(50) 


Prom  the  solutions  of  the  E\  equation,  as  described  by  (25)  through  (28), 
we  have  determined  that  C2  and  C3  are  functions  of  the  input  Co  and  the  perturbation 
term  b,  (c.f.  (14))  so  (46)  is  expressed 


D,  =  (yi'r 3  +  y2-rs)C,  + 


(51) 


J/o  Po 

We  now  make  the  observation  that  since  Di  is  the  value  of  Eq{x,  z)  at  x  =  0  and 
since  we  are  representing  the  total  wave  propagating  in  the  waveguide  by  (11),  then 
we  have 

D1  =  Ay  (52) 
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However,  we  already  have  expressed  b  in  terms  of  Ay  and  a  through  (14)  and  since  a 
has  been  found  independently  by  (37)  then  we  may  substitute  these  results  into  (51) 
and  (52)  to  determine  an  expression  for  Ay  in  terms  of  the  input  Co  as 


Ay  - 


(Vi  •  r3  +  y3  •  r5)C0 


Vo  -  (Vi  -U  +  ys-  r6)(^4je 


(53) 


Note  that  if  we  ignore  the  b  terms  in  (26)  and  (27)  since  these  will  be  much  smaller 
than  the  input  wave,  then  we  obtain 


A  _  (Vi  ’  r3  +  P2  •  rs) 

y  yo  0 

Then  the  efficiency  is  found  from  (31)  and  (45)  as 

0  AvAl  n  J,m{yi-r3  +  y2-ro)(yi-r3  +  y2-r5y 

r\  =  2  ad  —  ^  *  =  2  ad - ; - 

V2  CqC0  rj2  yoyo 


(54) 


(55) 


IV  Illustration  of  Results  and  Interpretation 

First  we  plot  the  variation  of  both  a  and  rj  as  a  function  of  thickness  t  for  a  given 
thickness  d  for  both  small  and  large  index  variation  in  Fig.  2a,  2b  and  2c  respectively. 
For  all  these  plots  a  gives  the  growth  of  the  field  (the  growth  of  the  power  would 
be  2a).  It  is  noted  that  the  efficiency  and  the  diffraction  do  not  track  each  other 
as  was  the  case  for  diffraction  out  of  the  guide  [8].  This  results  from  the  matching 
procedure  and  the  presence  of  standing  waves  in  the  structure.  The  result  for  77 
shows  twice  as  many  peaks  as  the  result  for  a  because  two  matching  procedures  are 
involved.  The  peaks  in  a  occur  for  t  =  A/2  and  the  peaks  in  77  occur  for  t  =  A/4.  In 
Fig.  2b  the  values  of  77  increase  dramatically  with  the  increase  of  index  change  as  one 
would  expect.  As  the  index  difference  increases  the  range  of  validity  of  the  calculation 
decreases  and  in  Fig.  2b  it  is  indicated  to  be  ~  0.07  fxm.  The  determination  of 
is  shown  in  Fig.  2c  where  we  plot  the  variation  of  the  effective  bounce  angle. 
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As  explained  previously  [8],  the  calculation  is  based  upon  the  constancy  of  (5q  the 
propagation  constant  in  the  waveguide  both  with  and  without  the  grating.  When  the 
effective  bounce  angle  goes  to  zero,  it  implies  that  the  loss  has  become  so  great  that 
no  effective  propagation  can  occur.  Therefore  in  the  case  of  a  large  index  change  of 
1.6  the  grating  should  be  etched  to  no  more  than  0.07  fxm  for  a  1.0  fim  guide.  In  Fig. 
3,  the  efficiency  is  compared  for  the  TEq  and  TEi  modes  for  the  weak  grating  case 
with  a  significant  apparent  increase  for  the  higher  order  mode  because  of  stronger 
grating  interaction. 

It  is  interesting  to  note  that  the  energy  is  being  diffracted  into  the  first 
order  wave  Eq  and  E2  (see  Fig.  1).  As  shown  by  (41)-(45)  the  second  order  wave  E\ 
in  the  grating  must  first  be  excited  to  provide  the  perturbative  driving  force  for  the 
Eq  and  E2  waves  and  the  wave  E\,  being  the  solution  of  (10),  is  excited  only  by  a 
grating  with  a  second  order  pitch.  This  means  that  Eq  and  _E2  would  not  be  produced 
if  the  grating  pitch  was  the  first  order. 

In  Fig.  4,  the  77  and  a  are  shown  as  a  function  of  the  tooth  width.  There 
is  a  broad  maximum  for  both  parameters  with  the  result  for  a  showing  symmetry  but 
the  result  for  77  showing  a  slight  displacement  to  values  of  w/A  <  0.5. 

In  Fig.  5,  77  and  a  are  shown  versus  guide  thickness  d  and  the  same 
behavior  is  noted  as  for  diffraction  out  of  the  guide.  However,  the  efficiency  shows 
more  oscillatory  behavior,  which  evidence  of  additional  standing  wave  effects  resulting 
from  the  second  matching  procedure  for  77. 
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V  Approximate  Second  Order  Grating  Pitch 


For  diffraction  into  the  guide,  we  need  to  consider  both  the  angle  9  and  the  grating 
pitch  A.  For  any  combination  of  parameters,  the  phase  matching  condition 

O7 T 

Po~-j  =  Pi  (56) 

must  be  satisfied  where  pi  is  the  component  of  the  input  wave  vector  in  the  2  direction, 
i.e. 

Pi  =  korii  sin  6 

Clearly,  for  non-zero  6 ,  then  P  ^  0  and  if  (56)  is  satisfied,  then  A  ^  As.0  where  As.0. 
is  the  second  order  value  for  the  chosen  wavelength.  Stated  alternately,  for  A  =  As  o., 
the  incident  wave  must  be  normal  to  the  guide  for  a  diffracted  wave  to  build  up  in 
the  guide.  Also  for  each  value  of  A  (0  As.0.),  there  is  an  unique  value  of  9  for  which 
diffraction  into  the  guide  will  occur,  which  increases  continuously  as  A  is  made  larger 
or  smaller  than  As  o..  Fig.  6  shows  the  situation  for  an  oblique  angle  9  and  a  grating 
pitch  A  7^  A s.o,.  Then  if  A  <  As.0.  the  phase  matching  condition  (56)  indicates  that 
the  diffracted  wave  propagates  in  the  opposite  direction  to  the  z  component  of  the 
input  and  if  A  >  As_0.,  the  propagation  is  in  the  same  direction  as  the  2:  component  of 
the  input  wave  (symmetrical  propagation  in  both  directions  can  only  occur  for  normal 
incidence).  These  comments  apply  to  approximate  second  order  grating  pitches,  and 
approximate  solutions  are  found  as  below.  For  third  and  higher  order  approximate 
grating  pitches,  there  will  be  more  than  one  value  of  9  to  satisfy  the  phase  match 
condition. 


Therefore  for  diffraction  into  the  guide  we  need  the  solution  to 
+  (fc2^2  _  01  )Ei  +  2iPi^  +  kl  BEq(z)  =  0 


This  equation  is  solved  as  in  the  out-of-the-guide  case  in  which  the  propagation  con- 


stant  in  the  grating  layer 


kg  ~  kgr  +  ikgi 

and  the  parameter  b  in  the  particular  solution  (given  by  (14)  for  ideal  second  order) 

b  bf 


now  have  real  and  imaginary  parts.  However  to  the  extent  that  a\  <C  Pi  we  can 
ignore  the  imaginary  part  so  that 

%  =  (kW  - 

b=  ,k°,„  n*BAse<“ 
k2  -  k'2 

The  remainder  of  the  solution  follows  the  calculations  of  (15)  through  (32)  to  produce 
modified  values  of  C\,  Cz,  Cz  and  Cq.  These  modified  values  are  used  to  predict 
modified  values  of  r2,  r8  and  a  by  comparison  with  (37). 

The  calculation  of  the  nonideal  second  order  grating  is  performed  in  the 
same  way,  with  the  alteration  in  (41)  -(45)  that  modified  values  of  C2  and  Cz  are 
used.  The  results  are  shown  in  Fig.  6  for  the  weak  grating  case  and  here,  for  the 
chosen  t  =0.2  fxm,  77  and  a  essentially  track  each  other. 


VI  Summary 

A  model  has  been  presented  to  describe  diffraction  into  a  waveguide  from  a  normally 
incident  beam.  A  closed  form  method  to  obtain  the  diffraction  efficiency  of  the  grating 
from  the  beam  to  the  guided  wave  has  been  described  that  can  be  implemented 
easily  on  a  PC  Platform.  The  calculation  shows  reversible  behavior  in  terms  of 
the  diffraction  constant  when  compared  to  diffraction  out  of  the  guide.  This  is  a 
result  of  the  conservation  of  energy  in  both  situations.  However,  the  efficiencies  are 
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dramatically  different  due  to  the  significant  energy  passing  through  the  guide  in  the 
present  situation. 

In  the  context  of  simple  phase  matching,  there  is  no  deviation  of  the  input 
angle  allowed  to  obtain  diffraction  i.e.  the  phase  matching  condition  states  that  if 
the  grating  pitch  and  incident  angle  are  not  matched  for  a  given  wavelength  then 
no  diffraction  into  the  guide  would  occur.  This  is  an  idealization  and  would  not  be 
expected  experimentally  because  of  parameter  variations,  specifically  the  refractive 
indices  which  are  generally  temperature  and  background  dependent.  It  is  noted  in 
our  discussion  that  the  particular  solution  for  the  wave  in  the  grating  region,  Egx, 
corresponds  to  the  calculation  of  the  coupling  coefficient  in  the  coupled  mode  theory. 


92 


References 


[1]  W.  Streifer,  D.  R.  Scifres,  and  R.  Burnham.  “Analysis  of  grating-coupled  radi¬ 
ation  in  GaAs:GaAlAs  lasers  and  waveguides”.  IEEE  J.  Quantum  Electronics, 
vol.  12,  no.  7,  pp.  422-428,  July  1976. 

[2]  K.  Ogawa,  W.  S.  C.  Chang,  B.  L.  Sopori,  and  Rosenbaum  F.  J.  “A  theoreti¬ 
cal  analysis  of  etched  grating  couplers  for  integrated  optics”.  IEEE  Journal  of 
Quantum  Electronics,  vol.  9,  No.  1,  pp.  29-42,  Jan.  1973. 

[3]  I.  I.  Mokhun’  and  M.  0.  Sopin.  “Geometric-optical  approach  to  the  problem  of 
excitation  of  a  corrugated  optical  waveguide  at  near-normal  angles  of  incidence” . 
Quantum  Electronics,  vol.  26,  No.  9,  pp.  836-838,  Sep.  1996. 

[4]  C.  Ghizoni,  B.  Chen,  and  C.  Tang.  “Theory  and  experiments  on  grating  couplers 
for  thin-film  waveguides”.  IEEE  Journal  of  Quantum  Electronics,  Vol.  12,  No. 
2,  pp.  69-73,  Feb.  1976. 

[5]  Y.and  Kubota  T.and  Takeda  M.  Nakaya,  T.and  Katoh.  “Diffraction  efficiency 
of  a  grating  coupler  for  an  array  illuminator  ”.  Applied  Optics ,  vol.  35,  No.  20, 
pp.  3891-3898,  Jul.  1996. 

[6]  J.  Y.  Duboz,  V.  Berger,  N.  Laurent,  D.  Adam,  and  J.  Nagle.  “Grating  cou¬ 
pled  infrared  modulator  at  normal  incidence  based  on  intersubband  transitions” . 
Appl.  Phys.  Lett.,  vol.  70,  No.  12,  pp.  1569-1571,  Mar.  1997. 

[7]  H.  Nishihara,  M.  Haruna,  and  T.  Suhara.  Optical  integrated  circuits.  McGraw- 
Hill,  1987. 

[8]  G.  W.  Taylor  and  C.  Kwan.  “  Determination  of  diffraction  efficiency  for  a  second- 
order  corrugated  waveguide  ” .  IEEE  J.  Quantum  Electronics,  vol.  33,  No.  2,  pp. 
176-186,  Feb.  1997. 


93 


[9]  G.  W.  Taylor  and  C.  Kwan.  “Correction  to  Determination  of  diffraction  efficiency 
for  a  second-order  corrugated  waveguide” .  IEEE  J.  Quantum  Electronics,  vol. 
33,  No.  6,  pp.  1041,  Jun.  1997. 

[10]  E.  G.  Loewen  and  E.  Popov.  Diffraction  gratings  and  application.  New  York:M. 
Dekker,  1997. 


94 


List  of  figures 


Fig.  1  Schematic  of  the  incident  wave  producing  transmitted  wave,  reflected  wave 
and  diffracted  wave  due  to  periodic  grating. 

Fig.  2a  Diffraction  efficiency  77  (dash  line)  and  diffraction  constant  a  (solid  line) 
versus  grating  thickness  for  a  rectangular  grating  (An  =  0.2). 

Fig.  2b  Diffraction  efficiency  77  (dash  line)  and  diffraction  constant  a  (solid  line) 
versus  grating  thickness  for  a  rectangular  grating  (An  =  1.6).  In  the  shaded 
region  t  >  tcr  and  the  calculation  is  not  valid. 

Fig.  2c  Diffraction  efficiency  versus  grating  thickness  for  a  rectangular  grating.  Solid 
line  present  TEq  mode  and  dash  line  is  TE\  mode. 

Fig.  3  Effective  bounce  angle  9eff  and  effective  waveguide  index  neff  versus  normal¬ 
ized  grating  depth,  t/d. 

Fig.  4  Diffraction  constant  and  diffraction  efficiency  versus  tooth  width  for  a  rect¬ 
angular  grating.  Solid  line  represents  diffraction  constant  a  and  dash  line  is 
diffraction  efficiency  77. 

Fig.  5  Diffraction  constant  and  diffraction  efficiency  versus  waveguide  thickness  for 
a  rectangular  grating.  Solid  line  is  diffraction  constant  a  and  dash  line  is  diffrac¬ 
tion  efficiency  77. 

Fig.  6  Diffraction  constant  and  diffraction  efficiency  versus  grating  pitch  for  the 
approximate  second-order  rectangular  grating.  Solid  line  is  diffraction  constant 
a  and  dash  line  is  diffraction  efficiency  77. 
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Appendix  A 

Particular  Solution  of  Wave  Equation 

The  x  dependent  differential  wave  equation  in  the  grating  is  obtained  from 
(12)  and  (13)  as 

rfiF  ~ 

+  k]Es  =  -kl  B  Aa(x)e-°  (Al) 

The  complementary  solution  (RHS=0)  is 

Egc(x)  =  C2eik9X  +  Cze~ik°x  (A2) 


the  particular  solution  is 

.2  5  Jk„x  f  A0(z)e~aze~ifc* 


—  OLZp—ikgT  ^  p 

0  -  c  **  /  - dx  -  kl  B  e~ik°x  /  "uv~/""r" — dx 

0  1  -2  ika  0  J  -2  ik0 


Egp  =  B  etK*x  J 


Ao(x)e  azetksx 


=  (-  2— 0  -)  B  Aye  az(cosKX  +  -sin/cr) 

K  Kg 


(A3) 


where  -2 ikg  is  the  Wronskian  and  A0(x)  in  (7)  has  been  used  for  the  grating  region 
to  represent  TE  modes.  Then  the  total  solution  is 

Egc(x)  =  C2eikaX  +  Cze~tksX  +  B  Aye~az( cos  kx  +  -  sin  kx)  (A4) 

y  k1  —  k“  K, 
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Appendix  B 


Solution  of  the  linear  equations  (15)-(24)  for  the  ideal  second  order  diffraction  condi¬ 
tion  yield  the  following  results  for  the  parameters  A,  B,C,  D. 


A  = 


8  8 

cos  (kgt)  +  —  sin  (kgt)  —  [cos(/ct)  +  -  sin(/ct)] 

K/g  K/ 


-^2-^)(1  +  “  *)  +  k9 £]  +  Yg  ~^  +  M]} 

"5(^)(1  ~~  jj^H008!2*2^  ~t)~  M  -  jr  sin[2;c2(^  - 1)  -  kgt}} 

-f  (cos  Kt  -\ —  sin  Kt){- - — )  cos [2^2 (d  —  <)] 

k  k>2  — 

+(j- cos  Kt  —  y~  sin Kt)(^2  —  f3)  sin[2Wd  —  i)l 

&2  k2  &2  —  ^3 


5  =  — sinffcot)  +  —  cos  (kgt)  —  [—  cos(Kt)  —  -7-  sin(Kt)l 
k2  y  kg  y  lk2  k2 

+  ^){sin[2fc2(d  -  <)  +  kgt]  -  y  cos[2k2(d  -t)  +  kgt ]} 
fcg^1  “  ^Hsin[2Mrf  -t)-  kgt ]  +  cos[2k2(d  -  t)  -  fc5t]} 
+(cos/c£  +  —  sinKt)(^2— ■  7 3 )  sin[2A;2(d  —  t)  ] 

AC  n^2 

—(•7-  cos  /ct  —  7-  sin/ct)("- — 7-)  cos^Wd  —  t)l 
&2  k2  k2-k3 


C  =  2(1  -  ^)cos(kgt)  -  (1  +  ^)(1  +  t£)(^2  +  ^3)cos[fcgt  +  2fc2(d  -  t)] 

n^2  »Cy  ^3 

-(1  -  ^)(i  -  cosCM  ~  _  0] 

P  =  2(^  ~  sin(fcgt)  -  (1 +  ^)(1  +  ^ )  sin[fcgt  +  2fc2(d  - 1)] 

+(1  ~  Fg^1  ~~  ^)(^T ^)sintM  -  2fc2(d-  *)] 
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Xq=  0.88  |im 
d  =  1.0  |um 
A  =  0.246  jim 
w/A  =  0.25 


Xq=  0.88  |im 


[g-OI.  x]u  ‘AON3IOIdd3  NOUOVaddlQ 
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Fig.  4 


DIFFRACTION  CONSTANT,  a  (cm-i) 


I 


